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1.0  INTRODUCTION 


1.1  FREE-STREAM  DIAGNOSTICS 

There  is  a  continuing  need  for  the  noninterfering  measurement  of 
both  static  gas  temperature  and  number  density  in  the  high-speed  flow 
fields  associated  with  aerospace  ground  test  facilities.  Specifically, 
there  are  immediate  applications  for  static  temperature  and  density 
measurements  in  free-stream  tunnel  flows,  for  the  studies  of  shock 
structure,  boundary  layers,  wakes,  turbulent  flow  fields,  condensation 
phenomena,  and  for  rocket  exhaust  plume  diagnostics. 

The  well-known  electron  beam  fluorescence  technique,  which  is 
capable  of  determining  these  parameters  (Ref.  1),  has  a  density  domain 
of  applicability  of  approximately  10^  to  10^  molecules/ cm^ .  The 
electron  beam  technique  is  density  limited  for  the  following  reasons: 

(1)  collisional  quenching  effects  on  the  measurements  increase  with 
density,  (2)  the  spatial  resolution  of  the  measurement  decreases  with 
increasing  density,  and  (3)  the  secondary  electron  effects  on  the  mea¬ 
surements  increase  with  increasing  density;  all  of  which  result  in  a 
decrease  in  accuracy  of  the  measurement  with  increasing  density. 
Therefore,  for  flow  fields  with  densities  greater  than  10^^  molecules/ cm^, 
static  gas  temperature  and  density  measurements  are  most  difficult. 
Laser-Raman  spectroscopy,  however,  provides  a  technique  for  static 
temperature  and  density  measurements  which  is  similar  to  the  electron 
beam  technique  in  that  the  measurement  is  local  and  nonperturbing  and 
possesses  molecular  specificity.  However,  Raman  spectroscopy  is  not 
susceptible  to  the  density  limitations  because  the  intermolecular  colli¬ 
sion  time  typically  exceeds  the  photon-molecule  scattering  interaction 
time  by  several  orders  of  magnitude. 

Complete  reviews  of  Raman  scattering  theory  may  be  found  in  Refs.  2 
through  5,  and  consequently,  only  a  brief  general  discussion  of  Raman 
scattering  will  be  undertaken  here.  Consider  an  incident  laser  beam  to 
consist  of  photons  of  energy  hrG.  On  collision  with  a  molecule,  the 
photons  may  be  elastically  scattered,  a  process  which  results  in  no 
change  in  the  molecular  internal  energy  states,  and  the  scattered  radi¬ 
ation  is  designated  as  Rayleigh  scattering.  Inelastic  collisions,  how¬ 
ever,  may  cause  the  molecule  to  undergo  a  quantum  transition  to  a  higher 
internal  mode  energy  level,  with  the  result  that  the  photon  loses  energy 
and  is  scattered  with  lower  frequency  (longer  wavelength).  The  encounter 
with  the  photon  may  also  cause  the  molecule  to  undergo  a  transition  to  a 
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lower  internal  mode  energy  level,  in  which  case  the  photon  is  scattered 
with  increased  frequency  (shorter  wavelength).  Figure  1  schematically 
depicts  the  transitions  occurring  in  a  molecule  for  Rayleigh  and  pure 
rotational  Raman  scattering,  and  Figs.  2  and  3  depict  a  typical  scat¬ 
tered  light  spectrum  that  might  be  observed  as  a  result  of  the  depicted 
transitions.  The  displacements  of  lower  frequency  have  historically 
been  referred  to  as  the  "Stokes  lines,  "  while  the  displacements  of 
higher  frequency  have  been  called  the  "anti- Stokes  lines.  " 
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Figure  1.  Partial  rotational  energy  level  diagram  depicting  transitions  resulting  in  the 
pure  rotational  Raman  spectra  and  the  Rayleigh  line. 
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Figure  2.  Calculated  pure  rotational  Raman  spectra  of  nitrogen  at  a 
temperature  of  300  K. 
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Figure  3.  Calculated  stokes  band  pure  rotational  Raman  spectra  of  nitrogen  and  oxygen 
in  air  at  a  temperature  of  50  K. 
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These  scattering  processes  may  be  visualized  as  the  perturbation 
of  the  molecular  eigenfunctions  by  the  electric  field  associated  with  the 
incident  photon,  inducing  transitions  via  intermediate  or  virtual  levels, 
and  the  process  occurs  with  characteristic  times  on  the  order  of  10"^^ 
to  10"^  sec.  It  is  the  shortness  of  this  characteristic  time  (rs)  which 
precludes  the  existence  of  gas  collisional  effects  at  reasonable  densi¬ 
ties;  i.  e.  ,  if  the  molecular  collision  time  (rc)  is  such  that  ts/tc  «  1, 
a  collision  is  unlikely  during  the  scattering  event,  and  no  collisional 
effects  are  to  be  expected.  For  a  rc  value  of  10“  ^  sec,  the  gas  density 
will  be  on  the  order  of  100  amagats,  and  for  tc  on  the  order  of  10 
sec,  the  density  will  be  on  the  order  of  10^  amagats.  Consequently,  no 
collisional  effects  are  to  be  expected  for  gas  density  values  less  than 
approximately  10^  to  10^  amagats.  This  is  to  be  contrasted  with  fluo¬ 
rescence  techniques  for  which  the  electronic  radiative  decay  processes 
exhibit  characteristic  times  (rr)  on  the  order  of  10"^  to  10  sec. 
Consequently,  collisional  effects  become  important  for  total  gas  den¬ 
sity  values  on  the  order  of  10"2  to  10"^  amagat,  or  lower,  depending 
on  the  intermolecular  collision  potential. 

For  diatomic  and  polyatomic  molecules,  the  relative  intensity  dis¬ 
tribution  of  the  pure  rotational  Raman  displacements  is  related  to  the 
rotational  temperature  which,  except  at  very  low  density,  is  equal  to 
the  translational  or  static  gas  temperature.  Therefore,  by  injecting  an 
intense,  monochromatic  laser  beam  into  a  gas  flow  field  and  spectro¬ 
scopically  analyzing  the  pure  rotational  Raman  spectrum,  the  rotational 
temperature  of  the  gas  may  be  determined.  Moreover,  this  measure¬ 
ment  (1)  will  be  independent  of  the  gas  flow  velocity  due  to  the  extremely 
short  scattering  time,  (2)  will  be  unaffected  by  gas  collision  processes, 
(3)  will  have  a  "point"  spatial  resolution  that  decreases  very  little  with 
increasing  density  and  requires  no  inversion  techniques  common  to 
absorption  measurements,  (4)  will  not  be  affected  by  secondary  particles 
since  there  will  be  no  ionization  process,  and  (5)  the  accuracy  will  in¬ 
crease  with  increasing  gas  density. 


1.2  INTENSITY  EQUATIONS 

The  observed  intensity  of  a  pure  rotational  Raman  line  of  a  diatomic 
or  linear  polyatomic  gas  is  known  to  be,  for  Stokes  scattering. 


RR 


n  virr  t(j) 


gj 


(2J  +  1) 


exp[-Bo.J(,J  +  l)hc/kTR] 


(1) 


12 


AEDC-TR-75-37 


where 


RR 


Y  = 
fi'  = 


K  = 


Tr(J)  = 

gj  = 
J  = 

%  = 

B0  = 

h  = 

C  - 

k  = 

Tn  = 


2 

the  rotational  Raman  cross  section,  cm  /ster 

_  g 

molecular  species  number  density,  cm" 
observed  scattering  volume 

the  solid  angle  subtended  by  the  collecting  optics 

the  laser  beam  intensity,  photons/ sec  cm^,  at  wavelength  (Aq) 

sensitivity  factor  of  optics,  spectrometer,  and  detector 

nuclear  spin  statistics  weighing  factor 

rotational  quantum  number 
molecular  rotational  partition  function 

rotational  constant  of  the  ground  vibrational  level 

Planck’s  constant 
speed  of  light 
Boltzmann’s  constant 
rotational  temperature 


The  J-dependent  portion  of  the  cross  section  may  be  separated  from 

aRRj;  L  e* ' 


RR 


CRRSJ 

A4 

AJ 


(2) 


where  the  line  strength  factor  (Sj)  is  given  by 


Sj  -  3(J  +  1)(J  +  2) /[ 2(2 J  +  1)(2J  +  3)] 


(3) 


the  wavelength  A  j  of  the  rotational  Raman  line  is  given  by 


Aj  =  l/[i7o  -  4Bo(J  +  3/2)] 


(4) 


and  vQ  is  the  laser  line  wavenumber  and  CpjR(cmb)  of  Eq.  (2)  is  a  con¬ 
stant. 
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If  one  defines  a  constant  C  as 

C  =  CRR  Ng  V  Q  Io 


(5) 


Eq.  (1)  may  be  written  as 


I 


J 


CTr(J)gj 


(J  +  1)(J  +  2) 
(2J  +  3) 


exp  [-  B0J(J 


+ 


l)hc/kTR] 


(6) 


Since  the  characteristic  rotational  temperatures  B0hc/k  for  the  molec¬ 
ular  species  investigated  in  this  study,  as  shown  in  Table  1,  are  much 
less  than  the  kinetic  temperatures  of  the  flow  fields  studied,  one  may 
make  the  following  approximation  for  the  rotational  partition  function: 


qr  =  TR  (aB0hc'k) 


(7) 


where  o  =  2  for  homonuclear  species  and  cr  =  1  otherwise.  This  approxi¬ 
mation  is  discussed  in  a  later  section. 


Table  1.  Molecular  Spectroscopic  Values 


Gas 

J 

— 

gj 

B0,  cm'^ 

B0hc/k,  K 

n2 

0,1, 2, 3,  ... 

1:  even 

2:  odd 

2.  001 

2.87 

02 

1,3, 5,7,  ... 

1 

1.438 

2.07 

CO 

0,1,  2,  3,  ... 

1 

1.  923 

2.  77 

co2 

0,2, 4,6,  ... 

1 

0. 3895 

0.  561 

From  Eqs.  (4)  and  (6),  we  find 


Pn 


Ij(kTR  aBohc) 

[Pp  -~4B  (J  +  3'2)]4  T  (J)g.ri  +  2) 

0  r  J[_2  (2J  +  3) 


=  Pn  C  - 


B  he 

O 

”k  Tn 


J(J  +  1) 


(8) 
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Written  more  conveniently  in  terms  of  relative  intensities. 


Equation  (9)  may  be  solved  by  iterative  computer  calculation  to  find  Tp 
using  measured  values  of  Ije^.  The  nitrogen  Tr  program  (PROGRAM 
RATROT)  for  the  CDC-1604B  computer  is  shown  in  Table  2. 


Table  2.  PROGRAM  RATROT 


TYPt  WtAt  NU * NUNuHm *  N^UMX Y , N SUMXSO 
TYPt  INTEGER  VU*vL 

OlMtNSION  I (24) *o(?4) *XY (?4) *XSQ(24) *DELTAY(24) *YT (?4) * 

•  X ( 24 ) *PtH(24) *WAvF (24) fTRR( 10*25)  * S 1 GMA T A ( 1 0 *  25 > * 

«  Tw(lOO),  Slb(lOo)'  S(100)*  SS(100)*  Y ( 1 o 0 ) * 

1  rY(l00)*  Z(lOO)*  w(loO)*  V»  Y  (100)  *  10(10)*  E  ( 3 1 )  *  0(31)  * 
]SZ(25)«  TPA  (25)  ♦ 

1  wtl  (40)  *f»dB  (40)  *Sh  (40) 

300  HEAL  333*  TwW,  LiLY 
IK  ( EOF  *  50 )  331  *33? 

331  SlOF 

332  LOLL-O 

00  334  L=l*?5 
r«(L)  =  T«W 
33**  CON  1  I NUE 
M  =  0 

1  Hfc  At  200  *RIJW 
LOLL*  LOLO  ♦  1 
N=  C 

M  =  M  ♦  1 

PHIM  100 
PHIM  103, HUN 
PHIM  104 
HtAl  20 1 *KMAX 
KKK  =  KM  A  X 
DO  3  K* 1  *  KMAX 
HEAL  202  *  I  (K ) 

IF ( (K*l ) /2.FQ.  K/2)  9999*3 
9999  I(K)=2«1(K) 

3  CONTINUE 
THU=200* 

4  sum  a  * o • 

SUMY  =  0 • 

SUMXY  =  (). 

SUMASUsO, 

SUMwsO. 

C  WAVELENGTH  CALCULATION 
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C  RfcLw  *  *2  ROT  AT iONAL  RAMaN  CALCULATION 

WO=E. 145E-05 
WfcjO* 1,/WO 
Hl*2 ,001 
UU  E  K=1 iKMflX 
Jj=*-1 

Wd  <K)s*HG-4.*HL* ( JJM.5) 

WAVt (K) = ( 1 ,/WH (K) )»l.tOR 
WBB(K)=(*BO/Wb(K))**3. 

SF(K>  =  (3,*(JJ*1)*(JJ^2) )/(?.* (2*JJ*3> > 

GIK)=SF  (K)/'fRG 

Y(K)=LOOF ( I (K)*Wbri(K)/G(K) ) 

X(K)  = 

X  Y (K)sK(K) * Y (K) 

XSO IK) =x (K ) *X (K) 

W(K)=I ( K ) ** 1 

SUMXsSUMX* (I ( K ) •* 1 ) *X ( K ) 

Sumy =sum  y ♦ <  r  <k>  **i >  *y  ik) 

SUMXY  =  SU*XY+ ( 1 (K) **  1  ) *XY (K) 

SUMXSG=5UPXSG*  < I (K) **i >**su (k) 

6  Sum »s sum w ♦ w ( k ) 

SUSl/MX  =  SUPX*SUMX 
RHOISUM=SLMX*SUMt 

SLOFE= (SUKW«SUMXr-RHOUSUM) / (SUMW*SUMXSQ-bUSUMX) 

OHDtEP= (SLMXSQ*SUMY-SUMX*SUMXY) / (SUMW*5UMXSQ-SQbUKX) 

rRs-2,87R/ SLOPE 

SUMUELSUsO, 

SUMmSU=0, 

SUMmXSQ=U. 

sum»sux=o. 
uo  e  k=i*kmax 
Y  I  (K)=SLOPE*X(K) ♦ORDCtP 
DfcLl AY (K)=Y (K)-Y f  (K) 

PtR(K)  =  lOO.*DELTAY (K) /YT  < K ) 

SOSLMW=bUPW*SUMW 

SUM»XSQ=S0M*X5U* ( I ( K ) **? ) «XSQ (K ) 

SUM»SUX  =  SUMWS(JX*  ( 1  (K)  **2>  *X  (K) 

summsossupwsq+i (k)**2 
fc  SUMCELSU  =  SUMUELSu^OEL  TAY (K) *deltay (K) 

SlGPAY=(SLMnELSO/ (KMAx-2«) )**.5 

S1GPaT= (?.879*SIumaY/ lSLOPE*SLOPE) > * ( (SCSUMW*SUMWXSG-2.*SUMW* 
*SUMX*SUM«SQX>SQSUMX*SUMWS(J) **.5)/ (SUMW*SUMXSQ-SUSLmX) 

TES1 TW=A4SF ( TR-TkG) /Tp 
IF  (TESTTP.LT. l,t-04)  8*7 

7  Thb=TR 
(30  10  4 

8  N=N*  1 

IF  (N.GT.l)  11*9 

9  DC  10  K  =  1,  KMAa 
JJ*K-1 

10  PH  IK  T  10S* JJ*I (K) *G(K) *Y(K) *X(K) *WAVE(K) *PER(K) 

PRINT  10b 

PRINT  107*PROL)SUm,N'SUmXY,SQSUMX,NSUMXSG 
PRINT  10H*SLOPE*URnCEP 
PRINT  109 

PRINT  110*KMAX,Th,SI(3MAT 
TRR(M’KMAX)  =  TP 
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S1GPATA(M,  KMAX)  a  SIuMAT 
KPAXaKMAA- l 
GO  TO  4 

11  PRIM  111 *KMAX*Th*SIGmAT 
TRR(M*KMAA)  a  TP 
SIGPATA(M*  KMaX)  a  SIGMAT 
KHAX=KHAX-1 

IF  (KMAX.GT.l)  4*400 
400  IF  (LULO.LT. LILY)  1*12 

12  CONTINUE 
KP  AX  =  KKK 
PRIM  65 

65  FORPAT  I1H) 

19  CONTINUE 

IF  <KMAX.C*T.1>  18*990 
18  CONIINUE 
PRIM  70 

00  13  J  =  1*  M 

PKlM  80*  TRR(J*kmaA),  TW(J>*  SlGMATA ( J*  KMAX) 

13  CONIINUE 

00  14  J  =  1*  M 
Y ( J )  =  TRR ( J*  KMaX)/Tw(J) 

S(J)  =  SIGMaTAU*  KMAX)  /  1 4  ( J) 

SS(*)  =  5 ( J) **2. 

YtL)  a  Y(J)/SS(J) 

Z(J)  =  l./SS<J) 

«(„)  a  SS<1)/SS(J) 

14  WY(w)  =  *(J)  *  *(J)  *  SS(J) 

SUMl  a  0.0  %  SUM2  a  O.o  S  SUm3  *  0.0  SSUm4  a  0.0 
1)0  15  J  =  1  ♦  M 
SUMl  a  SUM  ♦  YY  ( J) 

SUMZ  a  SUP?  ♦  Z(J> 

SUM3  a  SUP  3  ♦  WYlj) 

15  SUM4  a  SUP4  ♦  W(J) 

THBAR  =  SUM1/SUM2 

SIGBAW  »  <SUM3/(SUM4  *  SUM4 ) )  **  0«5 
PHIM  30*  TRBAR*  SIGBa« 

PRIM  210.  KMAX*  M 
SSS  a  SIG8AR  *  T mm) 

SZ (KMAX )  *  1 .  /  (SSS*SSS> 

TRA (KMAX )  a  TRBAW*T*W*SZ(KMAX) 

KP AX  a  KMAX  -  1 
PRIM  ?3 

23  FORMAT  (1R0*  ///i 
GO  TO  19 

990  CONTINUE 
PRIM  ?3 

IF^KKK.GT.a1  991*993 

991  CONIINUE 
PRIM  1000 
PRIM  1003 

DO  992  UM  a  4*KKK 

51  a  0.0 

52  a  o.O 

DO  994  J  *  3*JM 
SI  *  SI  ♦  TRA(J) 
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994  S2  =  S2  ♦  S/(J) 

THAtS  =  S1/S2 
S 1  GABS  =  ( 1 ,/S2) **Q.5 
SIGABS  =  ?.*SlGAbS 
C  ThL  2  SIGMA  ERROk  IS  uUOTeq 

992  PH  I M  1 00 1 ♦ JM» TRABS* S IGAbS 
GO  10  300 

993  PRIM  1002 
GO  10  300 

1000  F0HNAT(1H  . l0X*30HwElbHT£n  AVERAGE  T  CALCULATION//) 

1001  FURMTUh  IS»SX*F  10.2»5^’F10*2/) 

1  002  FUNf'flTUr  *3X.26FiFOUR  LINF  VALUF  IS  AVERAGE/) 

1003  ECR^ATdh  .gX.HHnO.L  InES.S*  ♦  1  0HAVE.  TEHR  .  ♦  5  A  ♦  9H  I  WC  SIGMA/) 

2in  FORMAT  (  l  F  0  ♦  5  A »  ^hUblNG.  1 3 »  )  8H  LINES.  AVERAGE  OF.  i2»  *SH  RUNS 
3u  FORMAT  ( 1 H  0 • 15HAVFRAGl  TE*P  =  .  Fl5»9*  5*«  1 6FtA  V£R  AGE  SIGMA  =  » 

1 F 10.5) 

7  u  FuRf'AT  (IF  .  ?  OHmF  a  SUhE  0  TEMPERATURE.  5X.  1 6Hw  ALL  TFMPERATURF. 

*I0X.  1  RRS  IAN  l)  ARD  UFVImIIUw  //) 

R  0  FORMAT  <1R  .  3*»  F  1  5  •  9  *  7  X ,  F15.9.  13X,  Flo. 5) 

100  FORMAT  (  IF  1 ♦/ .47a. PSHHOTAT IONAL  BOLTZMANN  PLOT) 

101  FORMAT  (IF  ./.57A.5HUblNG) 

102  FORMAT  (IF  .  /  ♦  52  A  .  1  4  Ft  M  UN  T  /  G-FACTOR) 

103  FORMAT  (If  ./  /  ♦  ?  1  A  .  ]2F)hUN  NUMBER  =.1X.A8) 

104  FORMAT  (IF  . // *  6A  ♦  l  RN  ♦  7X  ♦  4H  I  (K  )  .  1  1  X  .4HG  (  K  )  ♦  15X  »4H'T  (K  )  ♦  15X  .  AH  A  (K  )  . 

*IJX*10H*avelENGTh.9X»27fi(oELTA  Y(K))/YT(K)»  PERCENT) 

105  FORMAT  (IF  *  4 X  *  I  2 ♦ 6 X » 1 5 ♦ 8 X ♦ E 1 1 • 4 ♦ R X . E 1 1 . 4 » R X . E 1 1  •  4  * P X  .  E  1  1  •  4  ♦ 
*12X.£11.h) 

106  FORMAT  (IF  ♦ / / ♦ 5 A ♦ i Jh bUM ( X ) * SUM ( Y )  . 6X . 9F N* SUM ( X Y )  . R X  ♦ 
*1JHSUM(X)«SUM1(X)  .6X.  luhN*SUM  (X*X)  ) 

1 0  7  FORMAT  (IF  .5X.E  1  1  ,4.oA  .El  1 .4.  7X .E 1 1 .4 . 7X.E  1 1  .4) 

108  FOHMAT  (if  *//.2iX  .7FISLOPF  =  »  1  X  ,  E  1 1 . 4  *  1 G  X  »  1  \  H I N  I  ER  CE  P  T  =  »1X»EU.4 

109  FORMAT  (  IF  ♦//) 

11J  FORMAT  (IF  .15X. 12.6H-LINFS.5X.24HR0TAT1UNAL  TEMPERATURE  =. 

*E 1 1 .4 ♦ sx ♦ 26FST ANuARD  DEVIATION  OF  TR  =.E11.4) 

HI  FORMAT  (IF  .  15X.  I?,6H-LINES.5X»24F(  s, 

*E  1 1  •  4  ♦  5X  .  2  6F<  *  ♦  E  l  1  •  4  ) 

200  FURmaT(AH) 

201  FORMAT  (  12 ) 

202  FORMAT  (15) 

333  FORMAT  (F 1 0 .3.  12) 
enu  trot 


The  program  essentially  performs  a  least-squares  fit  of  the  data  points 
to  a  straight  line,  and  the  slope  of  that  line  is  used  to  calculate  the  ro¬ 
tational  temperature  Tr.  The  iteration  process  continues  to  calculate 
a  temperature  until  consecutive  calculated  values  are  within  0.01  per¬ 
cent.  It  should  be  noted  that,  in  the  least-squares  fit,  the  line  inten¬ 
sity  values  (Ije^)  are  weighed  according  to  the  reciprocal  of  the  square 
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of  the  fractional  standard  deviation  of  the  relative  line  intensity.  There¬ 
fore,  the  more  intense  and,  thereby,  more  precise  lines  of  the  spec¬ 
trum  are  directly  given  more  weight  in  the  temperature  determination 
process.  For  details  of  the  temperature  determination,  see  Appendix  A. 

From  Eq.  (1)  the  total  observed  intensity  of  the  pure  rotational 
Raman  Stokes  lines  is 


/photons\  _  V^'c 
\  sec  / 


-  j?0<TRRJTr(J>8.l(2J  +  «  exP  t  -B0J(J  +  l)hc'kT„] 


(10) 


To  avoid  the  measurement  of  absolute  intensities,  it  is  noted  that,  if 
li^e*  is  measured  for  a  known  static  gas  density  and  temperature,  a 
calibration  factor  (Cp)  can  be  determined  such  that  the  number  density 
(Ng)  of  a  gas  flow  can  be  determined  from 


N  = 


CF  Ip  gj(2j  +  1)  exp  [-B0J(J  +  l)hc/kTR] 

j?o  gJ  [■  '7IJT 3)~]  {exp  ["  BoJ(J  +  1)hc/kTR1} 


in  which 


(ID 


Jmax  =  rotational  quantum  number  corresponding  to  the  largest  wavelength 
change  Wiich  is  detectable 

Ijel  =  total  relative  intensity  of  the  pure  rotational  Raman  Stokes  lines 
Tpp  =  the  measured  rotational  temperature 

^max 

qr  2^  gj( 2 J  +  1)  exp  [ —  B0J(J  +  DhcTTp] 


CF 


m  ax 
2 

J=o 


8J 


['rnfa)  2]{»P  +  Dhe/kTjj- 


re 

T 


^max 

gj(2j+  1)  exp  [-B0J(J  +  l)hc/kTj 


(12) 
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In  Eq.  (12),  Nc  and  Tc  are  the  known  static  gas  density  and  tempera¬ 
ture,  respectively,  and  IrJ^  =  calibration  value  of  the  total  relative 

intensity  of  the  pure  rotational  Raman  Stokes'  lines.  The  calculations 
indicated  in  Eqs.  (11)  and  (12)  are  performed  by  PROGRAM  NUMRHO 
as  shown  in  Table  3. 


Table  3.  PROGRAM  NUMRHO 


*►  8HTRAN  LS/UO 


It 

DIMENSION  HI  <25)#KT(23>#KN<25> 

St 

HtAL  10*  NO 

j: 

RtAD( 105*7>n 

4t 

7 

►  OHMAH  12) 

5t 

111 

RtADUU5#l)  18*  TO#  NO 

6  t 

1 

FORMA) ( 3E 10*3) 

/  I 

WfcAD( 105*2) AMAX 

at 

2 

FONP1A1  1  12) 

*t 

Rt  AD ( 105# 3 ) ( K 1 (K)#NT< A)#A«1#KMAX> 

lot 

3 

KONMAM2E10*3) 

lit 

(,►  •  NO/  1  0 

let 

HA1p2*0/3*0 

ut 

Rd  1  •  1  •  0 

I4i 

DO  101  J- 1  *  23*  2 

1st 

KJ4  J 

1 6  t 

KA 1  PRA 1  ♦  (  (  K  J*  1  «  JMKJ42*  )/  l2.*KJO*  )  *fcXH  (  *2 . 8  79-R  J*  (  *  J*1  .  >/T8  )  )/2#Q 

in 

Kdl *Koi ♦ ( (2««KJ41« )*EAK(-2»879»RJ»(RJ4l» )/T0) )/2*0 

lot 

101 

CONT INOt 

i*t 

DO  JOS  U*2*  24*  2 

2ot 

HO-  J 

21 1 

KA1«MA14(MJ41« )MKJ*2* )/(2«»HJ43* )4EAH(.2#879»HJ«(KW41« )/!8) 

22 1 

HB1*ko14(2«*KJ41. )»EXH(»2*8;9*HJ*tHj4i« )/T0> 

23 1 

105 

CONTlNUt 

24  t 

<:►  -CF-HAl/Kdl 

2D  t 

WHITE t 108# 112 >NAl#Kdl# Of 

261 

lie 

F0HMA1 (1H1#»KA1  ■•# lPt 1 3 • 6# 5X# 8H01  ■ *# 1 Pt 1 3 • 6# 5X # *Cf  ••»lPtl3*6> 

2;t 

DO  102  A  •  1  *  UMAX 

28  t 

Hd-2*/3« 

2*1 

HA ■ 1*0 

30  t 

DO  103  J« 1 #  23*  2 

311 

HJ- J 

32  t 

KA«HA4(  (2.4KJ  +  1. )«ExHC2#8794Hj*(H04lt )/KT(K)  )  )/2.Q 

33  t 

ROpKOM ( HO*  J • )*(HD42. ) / l 2«  *H J*3 • )»£XHi-2«87*-HD*(KJ-l« )/Hf(K) ) )/2* 

34  t 

103 

CONTlNUt 

3d  i 

DO  104  0*2*24*2 

36  t 

RU'D 

3M 

HAPHA*12»*HJ41« )*tAPt-2*879*RJ* (KJ*1. )/KT  <<>) 

38  t 

Hd-Hd^tNJ*!, )*(RJ*2. )/l2.*HU43.  )-EXP(«2»879»HJ*{Hj41.  )/rTIA) ) 

3*1 

104 

CONTlNUt 

4o: 

KN( A ) «C*  *KI  (A  )*HA/Kd 

M : 

WHI  TL  1 1UB*  1 14  )A 

42» 

114 

►  OH'iA  H  1H0#  *A  p*#^) 

43» 

pK{Tt  Vl08*U31HA*Hd 

44  1 

113 

► OWWA U 1«0* *HA  ■•* 1PE13«6* 5X* *Htt  -**lPtl3»6) 

43  ; 

102 

CONTlNUt 
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46: 

WKJ  rt< 1U8/6) 10*  TB*N0 

4/: 

6 

FOKMAr(lHJi»lO  -»/ 1PL13«6*5X* *T»  1 HE  1 3 • 6 *  5 X* *N0 

48  t 

WKJTtl 108*5) 

4*1 

5 

F OKM  A | ( 1HU#»  10X,*i (*>*/ 14X#*T (*>*# 14X,*N{K)*) 
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Carbon  dioxide  (CO2),  although  a  polyatomic  species,  in  its  ground 
electronic  and  vibrational  state  may  be  treated  as  a  linear  molecule 
with  no  electronic  angular  momentum.  Unlike  N2,  O2,  and  CO,  at 
room  temperature  the  twofold  degenerate  deformation  vibrational  mode 
has  approximately  8  percent  of  the  CO2  molecules  in  the  first  vibra¬ 
tional  level,  and  consequently,  Eqs.  (4),  (6),  et  seq.  to  be  rigorously 
correct  should  explicitly  indicate  that  B0  and  are  required.  How¬ 
ever,  in  light  of  the  spectral  resolution  employed  in  this  study,  differ¬ 
ences  in  the  rotational  constant  Bv  due  to  rotation-vibration  coupling 
are  not  observed  and  will  henceforth  be  neglected. 

The  pure  rotational  Raman  spectrum  of  O2  is  more  complex  than 
that  of  N2  or  CO,  for  example.  The  ground  electronic  states  of  N2  and 
CO  are  both  i.  e.  ,  singlet  states,  but  the  ground  electronic  state 
of  O2  is  a  3T,  or  triplet  state.  This  difference  manifests  itself  as 
additional  structure  of  the  rotational  energy  levels  and  hence  of  the 
rotational  Raman  spectra  as  two  satellite  components  adjacent  to  pri¬ 
mary  Raman  lines.  As  will  be  seen  later,  however,  with  moderate 
spectral  resolution,  the  satellite  lines  are  not  resolved,  i.  e.  ,  the  in¬ 
tensities  of  the  satellite  lines  and  primary  line  are  approximately 
summed  together.  Since  the  sum  of  the  line  strengths  for  the  satellites 
and  primary  line  are  the  same  as  the  line  strength  for  the  rotational 
Raman  line  of  a  molecule  with  a  ground  electronic  state,  we  are, 
therefore,  justified  in  cases  of  moderate  spectral  resolution  in  treat¬ 
ing  O2  as  a  molecule  with  a  ground  electronic  state. 

Spectroscopic  constants  and  statistical  weights  for  use  in  Eqs.  (9) 
and  (11)  for  the  molecular  species  studied  are  given  in  Table  1. 
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2.0  EXPERIMENTAL  APPARATUS 


2.1  LABORATORY  SYSTEM 

For  laboratory  verification  studies,  the  experimental  arrangement 
shown  in  Fig.  4  was  used.  The  laser  source  was  a  Coherent  Radiation 
Model  52 -B  argon  ion  laser  nominally  operated  at  a  power  level  of  1.  5 
watts  at  either  a  488.  0-  or  514.  5-nm  wavelength.  As  shown  in  Fig.  4, 
the  laser  beam  was  deflected  90  deg  by  a  front  surface  mirror  and 
focused  by  a  50.  8-mm-diam  fused-silica  lens  of  250-mm  focal  length 
into  the  center  of  the  scattering  chamber.  By  using  the  relations  of 
Schwiesow  (Ref.  6),  which  are  given  in  Appendix  B,  the  diameter  of 
the  observed  scattering  volume  was  220.  5  Pm.  The  scattering  cell 
was  designed  similarly  to  that  of  Gerry  and  Rose  (Ref.  7),  and  it  was 
made  of  black  anodized  aluminum.  The  laser  beam  entered  the  scatter¬ 
ing  chamber  through  a  fused-silica  window,  and  from  the  focal  point 
within  the  chamber,  the  beam  expanded  and  entered  a  laser  dump. 

This  dump  consisted  of  an  offset  5 -deg  cone  of  72 -cm -length  black  ano¬ 
dized  aluminum  fitted  inside  a  closed  end  piece  of  black  anodized  alumi¬ 
num  tubing  which  was  bolted  to  the  scattering  cell.  The  focal  region  of 
the  laser  beam  was  observed  at  90  deg  through  a  fused-silica  window. 
Directly  across  the  cell  from  the  observation  port  was  a  viewing  dump 
that  consisted  of  a  black  anodized  15-deg  cone  of  25.4-cm  length  housed 
in  a  closed-end  piece  of  black  anodized  aluminum  tubing  which  was 
bolted  to  the  scattering  cell.  This  viewing  dump  provided  a  dark  back¬ 
ground  for  the  observation  of  the  laser  beam. 

The  entire  scattering  chamber  was  evacuable  to  10  torr  and  was 
essentially  leak-free.  Gases  of  interest  were  leaked  into  the  chamber 
through  a  Matheson  Model  G134  filter  and  through  a  series  of  valves 
including  a  needle  valve  which  were  connected  to  the  high-pressure  gas 
bottles.  Pressures  in  the  scattering  chamber  were  measured  with  a 
calibrated  Wallace  and  Tiernan  gage,  and  chamber  temperature  was 
measured  with  a  copper- constantan  thermocouple  connected  to  an 
L  &  N  temperature  indicator.  The  chamber  temperature  could  be 
lowered  by  flowing  LN2  through  insulated  copper  coils  wrapped  around 
the  scattering  cell. 
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The  focal  volume  of  the  laser  beam  was  61 -cm  from  the  7.  6 -cm 
diam,  30.  5 -cm  focal  length,  fused-silica  collection  lens.  Collected 
light  was  focused  onto  the  entrance  slit  of  the  spectroscopic  element 
located  61  cm  from  the  collection  lens.  A  magnification  of  one  was 
thus  achieved.  The  spectroscopic  element  was  a  Spex  1402  double 
spectrometer  with  a  Czerny-Turner  mount.  The  two  ganged  102-  by 
102 -mm  gratings  with  1200  grooves/mm  were  blazed  at  5000  A  with 
an  efficiency  of  approximately  70  percent,  and  silvered  mirrors  pro¬ 
vided  maximum  throughput.  The  focal  length  was  0.  85  meters,  and 
the  reciprocal  linear  dispersion  was  approximately  0.44  nm/mm. 

The  aperture  was  approximately  f/ 8  and  was,  therefore,  matched  to  the 
collection  optics.  The  spectrometer  could  be  either  synchronously 
scanned  at  rates  of  from  0.2  to  100  nm/min  or  digitally  (incrementally) 
scanned  at  rates  of  from  0.005  to  125  nm/min.  Radiation  passing 
through  the  exit  slit  was  collected  by  a  47.  6-mm-diam,  25.  4-mm  focal 
length  glass  lens  and  focused  onto  the  photocathode  of  the  photomulti¬ 
plier.  A  schematic  of  the  collection  and  spectrometer  optics  is  shown 
in  Fig.  5. 
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The  photomultiplier  tube  was  an  EMI-9502B  with  an  S-ll  spectral 
response  (see  Fig.  6)  and  a  10-mm-diam  photocathode,  and  the  tube  was 
contained  in  a  thermoelectrically  cooled  Products  for  Research,  Inc.  , 
Model  TE-104  housing.  The  housing  could  operate  continuously  at  -26°C, 
which  provided  a  reduction  in  dark  current  by  approximately  a  factor  of 
45  from  the  room  temperature  value.  The  schematic  of  the  Products 
for  Research,  Inc. ,  photomultiplier  voltage  divider  network  is  shown  in 
Fig.  7.  Power  to  the  photomultiplier  was  supplied  by  a  well-regulated, 
highly  stable  Fluke  405 B  power  supply. 


Figure  5.  Schematic  of  collection  and  spectrometer  optics. 
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Figure  6.  Spectral  response  curve  of  EMI-9502B  photomultiplier. 
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Unless  Otherwise 
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Resistors  220  K,  2%,  1  w 
Capacitors  0.01  jj.f ,  1  kv 
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Figure  7.  Voltage  divider  network  of  EMI-9502B  photomultiplier. 
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The  output  of  the  photomultiplier  tube  was  processed  by  an  ORTEC 
photon  counting  system  which  is  shown  in  the  block  diagram  of  Fig.  8. 
The  coaxial  output  cable  from  the  photomultiplier  tube  was  terminated 
with  a  50 -ohm  resistor  at  the  Model  454  amplifier  input.  The  ampli¬ 
fied  pulses  were  then  passed  through  a  Model  436  discriminator  to  both 
a  Model  715  dual  counter/timer  and  a  Model  441  ratemeter.  The  rate- 
meter  output  was  displayed  by  either  a  Clevite  Brush  Mark  250  or  a 
Honeywell  Electronik  19  strip- chart  recorder. 


Figure  8.  Block  diagram  of  photon  counting  system. 


2.2  RESEARCH  CHAMBER  SYSTEM 


For  the  flow  measurements  in  the  von  Karman  Facility  4-  by  10 -ft 
Research  Chamber,  the  experimental  arrangement  shown  in  Fig.  9  was 
used.  The  vacuum  chamber  was  a  cylindrical  stainless  steel  vessel 
4  ft  in  diameter  and  approximately  10  ft  long  and  is  shown  in  Fig.  10. 
Initial  pumping  capabilities  were  provided  by  an  8500  liter/min  mechani¬ 
cal  pump  for  rough  pumping  and  15 -cm  oil  diffusion  pump  for  inter¬ 
mediate  pumping.  Final  pumping  was  provided  by  both  77  K  liquid  Ng 
and  20  K  gaseous  He  liners  inside  the  chamber.  With  no  gas  flow,  an 
ultimate  pressure  of  10"^  torr  was  achieved,  and  with  mass  flow  rates 
exceeding  1  gm/sec,  background  pressures  less  than  2  x  10 "3  torr 
could  be  maintained.  Both  ionization  and  alphatron  gages  monitored  the 
chamber  pressure.  The  gas  sources  used  were  either  a  sonic  orifice  of 
exit  diameter  of  1.  325  mm  or  a  conical  nozzle  of  1.  04-mm  throat  diam¬ 
eter,  exit  plane  area  ratio  of  13.4,  and  half-angle  of  approximately 
14.  5  deg.  Schematics  of  these  sources  are  shown  in  Figs.  11  and  12. 
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Figure  10.  Schematic  of  4-  by  10-ft  Research  Chamber. 
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Figure  11.  Schematic  of  sonic  orifice. 


Conical  Nozzle 


Figure  12.  Schematic  of  conical  nozzle. 

The  sources  were  mounted  on  an  x-y-z  motor-driven  traversing  mech¬ 
anism  to  provide  flow  field  profile  studies  with  fixed  optical  instrumen¬ 
tation.  Imprecision  of  the  axial  movement  was  0.013  cm.  Reservoir 
temperatures  were  determined  with  copper -constantan  thermocouples, 
and  reservoir  pressures  were  obtained  with  either  an  absolute  pres¬ 
sure  transducer  or  gage.  Gases  were  supplied  to  the  flow  generators 
from  a  high-pressure  K -bottle  through  a  series  of  valves  and  a  needle 
valve  to  provide  fine  control  of  the  reservoir  pressure.  The  gases  and 
purity  grades  used  were  as  follows:  N2,  prepurified,  99.  998 -percent 
purity;  O2,  industrial  grade,  99.  95 -percent  purity;  CO,  C-P  grade, 

99.  5 -percent  purity;  CO2,  Coleman  grade,  99.  99-percent  purity.  To 
minimize  effects  of  dust  and  particulate  matter  from  the  bottles  and 
lines,  two  Millipore  filters  were  installed  in  parallel  upstream  of  the 
flow  generators.  These  filters  removed  particles  of  diameters  of 
0.  025  ijl m  and  larger.  A  schematic  of  the  gas  inbleed  system  is  shown 
in  Fig.  13. 
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Figure  13.  Schematic  of  4-  by  10- ft  Research  Chamber  gas  inbleed  system. 

Once  again,  the  argon  ion  laser  was  the  light  source.  The  normal 
vertical  polarization  of  the  laser  beam  was  rotated  90  deg  by  a  half¬ 
wave  plate,  and  the  1.4-mm-diam  laser  beam  was  directed  into  a  28X 
beam  expander.  The  resulting  3  9.  2-mm-diam  beam  was  focused  by  a 
50.  8-mm-diam,  1000-mm  focal  length,  fused-silica  lens  and  deflected 
upward  through  a  fused-silica  window  into  the  center  of  the  test  cham¬ 
ber  by  a  front  surface  mirror  as  shown  in  Fig.  9.  By  using  the  rela¬ 
tions  given  in  Appendix  B,  the  diameter  of  the  scattering  volume  was 
31,5  Mm.  As  the  laser  beam  defocused,  it  entered  a  laser  dump,  the 
entrance  of  which  was  a  black  cardboard  tube  which  extended  into  the 
chamber  to  within  22.  5  cm  of  the  chamber  centerline.  The  latter  stage 
of  the  dump  was  a  black  anodized  aluminum  tube  capped  with  a  black 
anodized  conical  section. 

The  focal  region  of  the  laser  was  located  at  the  focal  point  of  the 
collection  optics  system  which  was  oriented  at  90  deg  with  respect  to 
both  the  flow-field  axis  and  laser  beam.  A  76.  2-mm-diam,  152.4-mm 
focal  length  quartz  lens  located  inside  the  chamber  collected  and  colli¬ 
mated  the  scattered  radiation  and  transmitted  it  through  a  glass  port  to 
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a  101. 6-mm-diam,  304.8-mm  focal  length  glass  lens.  This  lens 
focused  the  scattering  volume  onto  the  entrance  slit  of  the  Spex  double 
spectrometer  with  a  magnification  of  two.  The  collection  optics  system 
external  to  the  chamber  was  shrouded  with  layers  of  green  felt.  For 
further  reduction  of  stray  light,  a  collection  optics  viewing  dump,  con¬ 
sisting  of  a  piece  of  closed-end  black  cardboard  tubing,  was  located  on 
the  side  of  the  test  chamber  collinear  with  the  scattering  volume  and 
the  center  of  the  collection  lens.  All  optically  accessible  surfaces 
inside  the  test  chamber  were  painted  with  a  flat  black  coating. 

The  signal  detection  and  processing  system  was  the  same  for  the 
flow  field  measurements  as  the  laboratory  studies. 


3.0  MEASUREMENTS 

3.1  LABORATORY  MEASUREMENTS 

Prior  to  laboratory  measurements,  a  careful  alignment  of  the 
scattering  volume  with  the  spectrometer  slit  had  to  be  accomplished. 
With  the  laser  operating  at  very  low  power,  a  visual  alignment  of  the 
laser  beam  with  the  scattering  chamber  centerline  was  completed,  the 
chamber  was  sealed,  and  noninhaled  cigarette  smoke  was  orally  blown 
into  the  chamber.  With  the  laser  now  operating  at  a  moderate  power 
level,  the  resulting  intense  scattered  light  allowed  easy  alignment  of 
the  image  of  the  scattering  volume  with  the  entrance  slit  of  the  spec¬ 
trometer.  Fine  adjustment  was  accomplished  by  using  the  Model  1450 
Periscope  Viewer  accessory  for  the  Spex  double  spectrometer.  With 
the  periscope  down,  the  image  of  the  scattering  volume  framed  by  the 
entrance  slit  was  projected  onto  a  viewing  screen.  By  lateral  adjust¬ 
ment  of  the  collection  optics  lens  and  slight  adjustment  of  the  spectrom¬ 
eter  legs,  the  scattering  volume  image  was  centered  precisely  in  the 
entrance  slit,  and  the  periscope  was  simply  retracted. 

_  c 

After  evacuation  of  the  chamber  to  10  torr,  the  gases  of  interest 
were  then  slowly  inbled  into  the  chamber  until  the  desired  pressure  was 
obtained.  The  Stokes  lines  of  the  pure  rotational  Raman  spectra  of  N2 
were  scanned  at  atmospheric  pressure  and  room  temperature  (see  Figs. 
14  and  16)  and  at  an  equivalent  room  temperature  pressure  of  10  torr 
and  temperatures  of  296,  215,  163,  and  104  K.  The  298  and  163  K 
scans  are  shown  in  Figs.  18  and  20,  respectively.  The  nitrogen 
Boltzmann  plots,  i.  e.  ,  the  plots  of  Eq.  (9),  are  shown  in  Figs.  15,  17, 
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19,  and  21.  The  measured  rotational  temperatures  determined  from 
slopes  of  the  straight  lines  are  all  within  ±2  percent  of  the  measured 
gas  temperature.  Scans  of  room  temperature  and  atmospheric  pres¬ 
sure  C>2  (Fig.  22),  CC>2  (Fig.  24),  and  CO  (Fig.  26)  were  made,  and 
the  respective  Boltzmann  plots  are  shown  in  Figs.  23,  25,  and  27. 
Again  the  rotational  temperatures  agreed  to  within  ±2  percent  of  the 
measured  gas  temperatures. 


Figure  14.  Pure  rotational  Raman  spectrum  of  nitrogen  at  atmospheric  pressure 
and  room  temperature. 


Figure  15.  Boltzmann  plot  of  pure  rotational  Raman  line  intensities  from  Fig.  14. 
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Figure  18.  Pure  rotational  Raman  spectrum  of  nitrogen  at  10-torr 
pressure  and  room  temperature. 


Figure  19.  Boltzmann  plot  of  pure  rotational  Raman  line  intensities 
from  Fig.  18. 
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N2  Rotational  Raman  Stokes  Band 
N  ~  3.38  x  1017  cm'3,  T  - 163  K 
200  pm  -  200  pm  -  200  pm  x  5  mm  Slits 
0. 1  nm/min  Scan  Rate 
~  200  counts/sec  Full  Scale 
Laser  Beam:  1. 9  w  at  514. 5  nm 


Figure  20.  Pure  rotational  Raman  spectrum  of  nitrogen  at  an  equivalent  room 
temperature  pressure  of  10  torr  and  a  temperature  of  163  K. 


Figure  21.  Boltzmann  plot  of  pure  rotational  Raman  line  intensities  from  Fig.  20. 
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Figure  24.  Pure  rotational  Raman  spectrum  of  carbon  dioxide  at 
atmospheric  pressure  and  room  temperature. 


Figure  25.  Boltzmann  plot  of  pure  rotational  Raman  line  intensities  from  Fig.  24. 
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Figure  26.  Pure  rotational  Raman  spectrum  of  carbon  monoxide  at 
atmospheric  pressure  and  room  temperature. 


Figure  27.  Boltzmann  plot  of  pure  rotational  Raman  line  intensities  from  Fig.  26. 
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It  became  obvious  that  minor  improvements  in  the  laser  and  col¬ 
lection  optics  could  decrease  the  minimum  number  density  from 
3.38  x  10^  cm_3  (equivalent  room  temperature,  10  torr)  to  nearly 
3.38  x  10*6  cm-^  (equivalent  room  temperature,  1  torr).  Further¬ 
more,  the  measurements  had  been  made  with  relative  ease,  and  suf¬ 
ficient  accuracy  was  attained  in  the  rotational  temperature  determina¬ 
tion;  therefore,  rotational  temperature  measurements  in  hypersonic 
flow  fields  were  initiated. 


3.2  FLOW-FIELD  MEASUREMENTS 

For  the  flow-field  measurements  in  the  4  by  10  Research  Chamber, 
not  only  was  a  careful  alignment  of  the  scattering  volume  with  the  spec¬ 
trometer  slit  required,  but  a  very  accurate  measurement  of  the  position 
of  the  scattering  volume  relative  to  the  free  jet  or  nozzle  exit  plane  was 
necessary.  With  the  laser  operating  at  moderate  power  and  focused  into 
the  approximate  center  of  the  chamber,  a  small,  flat-faced,  incandescent 
light  with  frosted  bulb  was  placed  opposite  the  scattering  volume  from  the 
collection  lens  with  the  front  surface  of  the  bulb  nearly  touching  the  finely 
focused  laser  beam.  Collection  optics  were  adjusted  in  such  a  manner 
that  the  image  of  the  lighted  bulb  was  in  focus  on  the  spectrometer 
entrance  slit  and  centered  on  the  scattering  volume  defining  aperture 
which  was  an  adjustable  iris  located  just  in  front  of  the  spectrometer 
entrance  slit.  With  two  strips  of  black  tape  placed  horizontally  across 
the  bulb  face  separating  the  desired  distance  for  good  spatial  resolution, 
the  defining  aperture  diameter  was  adjusted  to  match  this  restricted 
bulb  image.  After  removing  the  bulb,  the  laser  was  operated  at  high 
power,  and  the  periscope  viewer  attachment  was  again  used  to  precisely 
align  the  scattering  volume  image  in  the  spectrometer  entrance  slit. 

A  K  &  E  telescope  mounted  outside  the  chamber  and  collinear  with 
the  collection  optics  was  focused  onto  the  scattering  volume.  After  par¬ 
tial  evacuation  of  the  chamber,  the  telescope  crosshairs  were  centered 
on  the  scattering  volume.  Axial  location  of  the  flow  generator  exit  plane 
was  then  readily  determined  by  axially  traversing  the  generator  until  the 
edge  of  the  exit  plane  and  the  telescope  vertical  crosshair  were  aligned. 
The  vertical  and  radial  zero  positions  were  measured  by  establishing  a 
flow  field,  making  vertical  and  radial  Rayleigh  scattering  intensity  pro¬ 
file  measurements,  and  determining  the  zero  positions  from  the  flow 
field  symmetry. 
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Initial  measurements  were  performed  for  sonic  orifice  expansions 
for  which  the  source  diameter  (D)  was  1.325  mm  and  the  orifice  wall 
thickness  (t)  was  such  that  D/t  >  25.  Figures  28  and  29  show  typical 
low  temperature  N2  rotational  Raman  spectra.  Spectrometer  scan 
rates  were  typically  0.  1  to  0.  2  nm/min  with  10-  and  5 -sec  time  con¬ 
stants,  respectively,  and  entrance,  middle,  and  exit  slit  widths  were 
200-200-200  pm,  respectively. 

N2  Rotational  Raman  Stokes  Band 

Sonic  Orifice,  D  =  1. 325  mm,  p0  =  3. 72  atm,  T0  =  290  K 

200 Mm  -  200 Mm  -  200 Mm  x  5  mm  Slits 

0. 1  nm/min  Scan  Rate 

~  1000  counts/sec  Full  Scale 

Laser  Beam:  1. 5  w  at  519. 5  nm 


Figure  28.  Typical  low  temperature  pure  rotational  Raman  spectrum  of 
nitrogen  in  a  hypersonic  flow. 


J(J  +  l) 

Figure  29.  Boltzmann  plot  of  pure  rotational  Raman  line  intensities  from  Fig.  28. 
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The  measured  axial  variation  of  the  temperatures  for  N2  stagnation 
pressures  of  2.  8,  3.  72,  5.  58,  and  7.  44  atm  are  shown  in  Figs.  30 
through  33,  respectively. 


r  =  1.4 

■Sherman-Ashkenas  Theory 
p0  =  2.79  atm,  T0  ■  287  K 
D  =  1.325  mm 
Prepure  N2 


Figure  30.  Axial  variation  of  TR/T0  for  a  free-jet  nitrogen  flow 
with  P„  =  2.79  atm. 


Figure  31.  Axial  variation  of  TR/TC  for  a  free-jet  nitrogen  flow 
with  =  3.72  atm. 
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Figure  32.  Axial  variation  of  TR  /T0  for  a  free-jet  nitrogen  flow 
with  PQ  =  5.58  atm. 
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Figure  33.  Axial  variation  of  TR/T0  for  a  free-jet  nitrogen  flow 
with  PQ  =  7.44  atm. 
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The  height  of  the  error  bars  indicate  the  2 -sigma  error  calculated  by 
PROGRAM  RATROT,  and  the  width  of  the  error  bars  indicate  possible 
ipositional  error.  For  X/D  values  less  than  9,  there  is  excellent  agree¬ 
ment  between  the  measured  values  and  values  calculated  using  the 
Ashkenas -Sherman  equations  (Ref.  8).  For  PQ  values  of  5.  58  and  7.44 
atm,  the  measured  rotational  temperatures  show  a  significant  increase 
over  the  calculated  values  for  X/D  greater  than  9.  For  these  high  pres¬ 
sures  and  axial  positions,  Rayleigh  scattering  measurements  have 
shown  massive  condensation  taking  place  (Refs.  9  and  10).  The  effects 
of  condensation  will  be  discussed  in  Section  4.  0.  The  lowest  number 
density  at  which  measurements  were  made  was  4.  52  x  10 cm'l 

Figures  34  and  35  show  typical  CO  and  O2  rotational  Raman  spectra 
taken  in  the  sonic  orifice  flow  field.  The  spectrometer  scan  rate  for  the 
CO  spectra  was  0.  1  nm/min  with  a  10 -sec  time  constant,  and  slit  widths 
were  150-200-150  pm.  The  axial  variation  of  the  CO  rotational  tempera¬ 
tures  for  PQ  =  3.  72  atm  is  shown  in  Fig.  36.  The  scan  rate  for  the  O2 
spectra  was  0.  2  nm/min  with  a  5-sec  time  constant,  and  slit  widths  were 
200-200-200  pm.  Axial  variation  of  O2  rotational  temperature  is  shown 
in  Fig.  .3  7.  Again,  there  is  excellent  agreement  between  measured  and 
calculated  temperatures  for  both  CO  and  O2  for  X/D  <  7.  However,  for 
larger  values  of  X/D,  the  measured  rotational  temperatures  show  a 
significant  increase  over  the  calculated  values,  and  this  effect  also 
correlates  with  massive  condensation  as  observed  by  Rayleigh  scattering 
(Ref.  10)  and  discussed  briefly  in  Section  4.  0.  The  lowest  number  den¬ 
sity  at  which  temperature  measurements  were  made  was  4.  54  x  10^  cm' 
in  O2  and  1.  34  x  10  ^  cm'^  in  CO. 

Figures  38  and  39  show  typical  low  temperature  CO2  rotational 
Raman  spectra  taken  in  the  sonic  flow.  The  spectrometer  scan  rate 
was  0.05  nm/min  with  a  5-sec  time  constant,  and  slit  widths  were 
50-100-50  Pm.  A  limited  axial  variation  (X/D  <  4)  of  the  CO2  rotational 
temperatures  with  PQ  =0.7  atm  is  shown  in  Fig.  40.  There  is  no  agree¬ 
ment  with  the  predicted  theoretical  variation  in  this  axial  region  for 
either  y  =  7/5  or  9/7.  Agreement  was  obtained  with  the  electron  beam 
results  of  Beylich  (Ref.  11)  using  appropriate  reservoir  scaling  laws 
(Ref.  12).  The  lowest  number  density  at  which  CO2  temperature  mea¬ 
surements  were  made  was  7.  97  x  10^6  cm'l 
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Figure  34.  Typical  low  temperature  pure  rotational 
Raman  spectrum  of  carbon  monoxide  in 
a  hypersonic  flow. 


Figure  35.  Typical  low  temperature  pure  rotational 
Raman  spectrum  of  oxygen  in  a  hyper¬ 
sonic  flow. 
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Figure  38.  Typical  low  temperature  pure  rotational 
Raman  spectrum  of  carbon  dioxide  in  a 
hypersonic  flow. 


Figure  39.  Boltzmann  plot  of  pure  rotational  Raman 
line  intensities  from  Fig.  38. 
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Figure  40.  Axial  variation  of  Tr/Td  for  a  free-jet  carbon 
dioxide  flow  with  PQ  =  0.70  atm. 

Measurements  of  axial  variations  of  both  rotational  temperature 
and  molecular  number  density  were  made  in  the  N2  flow  field  of  the 
conical  nozzle  for  PQ  values  of  3.41,  6.82,  and  10.2  atm.  The  results 
are  shown  in  Figs.  41  through  46.  The  number  density  measurements 
are  in  excellent  agreement  with  the  MOC  calculation  of  the  flow -field 
density  over  the  entire  X^/Dt  range  of  measurements.  For  the 
PQ  =  3.41  atm  condition,  the  measured  rotational  temperatures  are  in 
good  agreement  with  the  MOC  calculated  temperatures  over  the  entire 
X^/Dt  range  of  measurements.  Significant  increases  of  measured 
rotational  temperature  over  the  calculated  values  are  observed  in 
regions  of  massive  condensation  as  determined  from  Rayleigh  scatter¬ 
ing  (Refs.  10  and  12)  for  the  PQ  =  6.  82  and  10.  2  atm  conditions.  The 
lowest  density  at  which  temperatures  were  measured  was  2.31  x  10^ 
cm~3. 

Measurements  of  radial  variations  of  both  rotational  temperature 
and  molecular  number  density  were  made  at  axial  positions  of 
X^/Df;  =  12.47  and  19.  79  for  the  PQ  =  10.  2  atm  condition.  Results  are 
shown  in  Figs.  47  through  50.  Number  density  measurements  are  in 
good  agreement  with  the  calculated  radial  profiles.  However,  the 
rotational  temperature  radial  profiles  are  in  significant  disagreement 
with  the  calculations.  Even  though  condensation  may  be  only  moderate 
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on  the  centerline,  Rayleigh  scattering  measurements  (Ref.  10)  have 
shown  more  condensation  off  the  centerline  at  these  axial  positions, 
and  this  may  well  explain  the  higher -than- calculated  temperatures  at 
positions  off  the  centerline. 


Figure  41.  Axial  variation  of  TH/T0  for  a  conical  nozzle 
nitrogen  flow  with  PQ  =  3.40  atm. 
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Figure  44.  Axial  variation  of  p/p„  for  a  conical  nozzle 
nitrogen  flow  with  Pc  =  6.80  atm. 
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Figure  45.  Axial  variation  of  TR/TG  for  a  conical  nozzle  nitrogen  flow 
with  Pe  =  10.2  atm. 
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46.  Axial  variation  of  p/pQ  for  a  conical  nozzle 
nitrogen  flow  with  PQ  =  10.2  atm. 
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.  Radial  variation  of  TR/TQ  at  Xt/Dt  *  19.79  for  a  conical 
nozzle  nitrogen  flow  with  P0  =  10.2  atm. 
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Figure  48.  Radial  variation  of  p/pD  at  Xt/Dt  *  19.79  for  a  conical 


Figure  49.  Radial  variation  of  TR/T0  at  Xt/Dt  =  12.47  for  a  conical 
nozzle  nitrogen  flow  with  PQ  =  10.2  atm. 


Figure  50.  Radial  variation  of  plpQ  at  Xt/Dt  =  12.47  for  a  conical 
nozzle  nitrogen  flow  with  PQ  =  10.2  atm. 
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4.0  ANALYSIS  AND  DISCUSSION 


4.1  MOLECULAR  PARTITION  FUNCTIONS 

Preliminary  to  the  analysis  and  discussion  of  the  results  for  Raman 
diagnostics  of  the  flow  fields  studied,  the  previously  mentioned  approxi¬ 
mation  of  replacing  qr  with  a  function  linear  in  Tr  will  be  discussed. 

The  approximation  of  Eq.  (7)  is  actually  twofold:  the  replacement  of  the 
summation  over  the  discrete  rotational  quantum  number  J  by  kTj^/  (crB0hc) 
and,  secondly,  neglect  of  nuclear  spin  degeneracy  factors.  The  first  of 
these  two  is  relevant  to  both  homonuclear  and  heteronu clear  diatomic 
molecules,  while  the  second  is  important  only  for  homonuclear  species, 
assuming  a  nonreacting  flow  field.  For  completeness,  the  spin  factors 
will  be  included  for  both  types  of  diatomic  species. 

First,  consider  a  heteronuclear  specie  AB  for  which  the  atomic  con¬ 
stituents  A  and  B  have  nuclear  spins  1^  and  Ig,  respectively.  Since  the 
spin  degeneracy  factor  is  (21^  +  l)(2Ig  +  1),  the  quantum  mechanical  par¬ 
tition  function  q^^  is  given  by 

q?M  =  (2Ia  +  1)(2I„  +  1)  2j  gj  expl-J(J  +  1WH  Tr]  /, 


where 


6  =  BQhc  /k 

is  defined  as  the  characteristic  rotational  temperature 
and 


gj  »  2j  +  1 


Since  B0  varies  as  the  reciprocal  of  the  moment  of  inertia  of  the  specie, 
increases  as  the  specie  reduced  mass  decreases.  The  classical  limit 
Cl  of  qQM  iS  given  by  the  limit  at  Tr/0j{  -►  ®,  or 

<£’  -  (21 A  +  1X21 „  +  1)  X  Tr  '<?r 
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It  is  to  be  noted  that  molecular  rotation-vibration  interactions  have 
been  neglected  so  as  to  allow  writing  B0  for  the  molecular  rotational 
constant.  This  neglect  is  not  serious.  Since 


N/Ng  =  (21 A  +  1)(2Ib  +  1)  x  gjexp[-J(J  + 
the  classical  limit  is  given  by 


(15) 


1  im 

T  R/ 0  p-»oo 


(Nj'Ng) 


gj  exP[-J(J  +  DVTr^Tr^r) 


(16) 


Therefore,  as  expected,  the  nuclear  spin  factors  are  eliminated  for  a 
nonreacting  sample,  and  one  can  effectively  use  Tr/0r  as  the  partition 
function  for  relative  number  density  calculations,  as  was  done  in 
Eq.  (7). 


For  homonuclear  species,  the  same  considerations  will  apply  inso¬ 
far  as  replacing  the  summation  of  q^^  by  a  simple  function  of  Tr,  but 
nuclear  spin  statistics  assumes  a  much  more  important  and  visible  role. 
Since  the  high  temperature  ratio  of  ortho-  to  para- molecular  species  is 
assumed  to  be  retained  during  the  expansion  process,  one  must  regard 
the  gas  as  a  mixture  of  two  inconvertible  species.  The  partition  func¬ 
tions  for  the  ortho-specie  q^^(o)  and  the  para-specie  q^-^(p)  are  given 

by 

q?M(o)  =  (21 A  +  1)(IA  +  1)  S  '  gj  exp  [- J(J  +  1)0r/Th]  (17) 


and 


q?M(p)  =  (21 A  +  l)(IA)S*gj  exp[-J(J  +  1)0r'Th]  (18) 

where  E'  designates  summing  over  (even  J;  odd  J)  integers  if  (-l)aP 
is  C+T;-l)  and  a  is  the  nucleon  number  of  the  atomic  constituent 
and  P,  the  electronic  parity  eigenvalue,  is  ±1.  Then,  £"  denotes 
summing  over  the  remaining  integer  set.  Since  the  following  are  true. 


lim  E  gj  exp  [  -  J(J  +  1)^r/Tr]  = 

TR  °ir°° 


lim  S 

tr/0r^° 


gj  exp  [  —  J( J  +  1)<9 H  Tr]  =  Tr  ;2 dR 


(19) 
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the  high  temperature  equilibrium  partition  function  q^M(E)  is 


qPM(E)  =  j(2I A  +  l)2Sj  gj  exp[-J(J  +  IWr'TJI/2 

and 

lim  q?M(o)  =  (21 A  +  1)(IA  +  1)(Tr/2  0r) 

1  R  ^  R^°° 


(20) 


(21) 


lim  qpM(p)  =  <21 A  +  1)Ia(Tr  2(9r) 

TR  0R^° 


(22) 


one  finds  the  high  temperature  ortho-  and  para-number  density  distri¬ 
butions  to  be 


lim  [Nj(o)'jg  -  29r|1^1)  gjexp(-J(j  +  WB  T„]/TR 

iR/aR->oo  \  A 


(23) 


and 


Tn  dR^°° 


iim  [Nj(P)/Ng]  =  2enl2i7njejexP[-JU  +  ,wr'tr]/tr  (24) 


respectively.  Obviously,  the  factor  20j{/{2lj\  +  1)  of  the  preceding 
equations  is  of  no  consequence  for  calculation  of  the  relative  number 
density  or  intensity  distribution,  but  the  spin-dependent  factors  do 
influence  the  observed  scattered  intensity  distribution.  For  N2, 

(-l)aP  is  +1  so  even  J  molecules  are  ortho  species  and  odd  J  are  para. 
IA  is  1  so  even  J  states  are  twice  as  prevalent  as  odd  J  species,  and 
the  observed  intensity  data  for  odd  J  scattering  have  been  multiplied  by 
(Ia  +  1)/Ia  =  2  to  allow  for  spin  statistics.  For  O2,  since  Ia  is  zero 
and  P  is  -1,  only  ortho  species  exist  and  scattering  from  even  J  is 
absent.  Similar  considerations  apply  for  CO2  as  for  O2. 

After  having  determined  the  high  temperature  limit  of  the  number 
density  distribution,  it  now  remains  only  to  define  what  constitutes  for 
practical  purposes  a  high-temperature  limit. 


Figure  51  shows  the  temperature  variation  of  the  ratio  QeO  of  the 
summation  of  gj  exp  [-J(J  +  1)0r/Tj{]  over  even  and  odd  J  values; 
obviously,  the  approach  of  this  ratio  to  within  the  experimental  precision 
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value  of  unity  defines  the  high-temperature  limit  of  the  particular  mea¬ 
surement.  It  is  seen  that,  for  Tr/  ~  4,  QeO  deviates  from  unity  by 
less  than  1  part  per  thousand.  Therefore,  by  defining  the  high- 
temperature  limit  as  Tr/Or  =  4,  the  molecular  constants  of  Table  1 
show  that  the  approximation  of  Eq.  (7)  is  valid  for  N2,  C>2,  CO2,  and 
CO,  so  long  as  Tr  >  11.4,  8.3,  2.2,  and  11.  1  K,  respectively.  These 
restrictions  are  satisfied  for  the  present  study.  It  should  be  noted  that 
for  H2>  although  not  included  in  this  work,  =  85.  4  K  so  Tr  must  be 
greater  than  341.  6  K  to  use  the  classical  value  of  qr.  Since  isotopic 
species  AalAa2,  ai  #  a2,  are  heteronuclear,  the  ortho -para  distinction 
disappears  and  the  even-odd  J  ratio  assumed  previously  is  no  longer 
valid.  However,  since  both  the  0^:0^  ratio  and  the  N^:N^  ratio 
exceed  100:1  no  isotopic  specie  correction  is  applied  to  the  previously 
derived  intensity  equations. 


Figure  51.  Temperature  variation  of  the  ratio  Qeo- 
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4.2  CONDENSATION,  RELAXATION,  AND  THEIR  EFFECTS 

Excluding  at  present  vibrational  relaxation  processes,  the  only 
significant  source  of  free -stream  anisotropy  is  that  of  condensation. 
Obviously,  for  a  gas  expanding  from  a  source  of  unsaturated  reservoir 
parameters,  some  region  of  the  flow  field  must  be  unsaturated  and, 
hence,  uncondensed.  The  task  at  hand  is  to  determine  the  spatial 
regions  of  the  flow  field  for  which  condensation  is  not  possible  and  for 
which  the  calculated  flow  parameters  should  agree  with  the  measured 
values.  For  an  isentropic  expansion,  it  is  known  that  the  static  pres¬ 
sure  P  varies  as 


y/ (y —  1 ) 

P  =  P0(T7T0)  Y  (25) 

and  7  is  assumed  constant.  Consequently,  as  Fig.  52  shows,  the  ex¬ 
pansion  definitely  proceeds  isentropic  ally  in  the  P-T  plane  until  inter¬ 
section  of  the  vapor  pressure  curve.  The  saturation  temperature  Ts 
is  then  used,  in  the  form  of  Ts/T0,  with  the  results  of  either  Ashkenas- 
Sherman  or  the  MOCS  to  determine  the  axial  location  of  saturation 
(X/D)s,  prior  to  which  point  the  flow  field  is  definitely  uncondensed  and 
beyond  which  point  the  flow  field  is  supersaturated  but  not  necessarily 
condensed.  The  vapor  pressure -temperature  data  of  Hilsenrath,  et  al. 
(Ref.  13),  were  used  for  all  species  studied  except  CO2  for  which  the 
data  of  Honig  and  Hook  (Ref.  14)  were  used.  The  saturation  locations 
(X/D)s  so  obtained  are  shown  on  the  appropriate  figures  displaying  the 
experimental  results  of  this  work. 


Figure  52.  Diagram  of  expansion  process  in  P-T  plane. 
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Since  the  condensation  process  itself  is  not  the  subject  of  this  work 
and  is  discussed  elsewhere  (Refs.  9  and  12),  only  the  manifestation  of 
the  phenomenon  as  it  influences  this  work  will  be  considered.  The 
axial  profiles  of  the  Tr/T0  results  for  N2  sonic  orifice  expansions,  as 
shown  in  Figs.  30  through  33,  exhibit  varying  degrees  of  agreement 
with  the  theoretical  predictions.  It  is  seen  that,  for  X/D  <  (X/D)s,  the 
agreement  of  the  data  with  the  prediction  of  Ref.  8  is  excellent.  For 
X/D  >  (X/D)s,  the  agreement  with  theory  of  the  PG  =  2.  7  atm  expansion 
results  of  Fig.  30  is  good,  but,  as  Figs.  31  through  33  show,  the  devi¬ 
ation  of  measured  and  predicted  values  increases  as  PQ  increases,  and 
this  increase  in  the  measured  values  of  Tr/T0  obviously  results  from 
the  heat  release  to  the  gas  due  to  the  exothermic  combination  of  mole¬ 
cules  in  the  condensation  process,  and  the  magnitude  of  this  increase 
is  dependent  on  condensate  mass  fraction.  Previous  Rayleigh  scatter¬ 
ing  results  (Ref.  9)  have  shown  that  condensation  onset  and  cluster 
growth  have  indeed  occurred  for  the  flow  fields  of  Figs.  30  through  33 
prior  to  detection  of  increases  in  Tr.  This  is  not  unexpected  since  it 
has  been  found  that  the  Rayleigh  scattering  measurements  are,  in 
general,  at  least  an  order  of  magnitude  more  sensitive  with  respect  to 
condensate  mass  fraction  than  is  the  rotational  temperature  measure¬ 
ment.  Identical  considerations  apply  to  the  O2  and  CO  results  of 
Figs.  36  and  37. 

Although  the  sonic  orifice  CO2  results  of  Fig.  40  have  been  pre¬ 
sented  previously  (Ref.  12),  the  Tr  data  are  included  here  for  com¬ 
pleteness  and  additional  analysis.  From  Fig.  40,  it  is  seen  that  Tr 
of  CO2  was  measured  over  the  axial  distance  range  0.38  <  X/D  <4.6. 
The  agreement  with  the  electron  beam  results  of  Beylich  is  to  be 
noted. 

Also  shown  in  Fig.  40  are  the  predictions  of  Ref.  8  for  the  self¬ 
similar  axial  dependence  of  expansions  for  gases  of  7  =  7/5  and  9/7. 

It  is  noted  that  the  axial  region  of  validity  of  the  predictions  for 
7  =  7/5  and  9/7  is  X/D  >  1.0  and  4.0,  respectively.  Obviously,  our 
data  show  no  agreement  with  these  predictions.  To  determine  if  con¬ 
densation  exists  for  this  particular  sonic  orifice  expansion,  we  have 
used  unpublished  isentropic  CO2  density  data  using  Rayleigh  scattering 
for  both  3.  0-  and  1.  325-mm  orifice  diameter  expansions,  and  Fig.  53 
shows  these  results  when  plotted  in  the  N-T  plane.  It  is  seen  that 
7  =  1.44  fits  the  isentropic  portion  of  the  expansion  and  that  deviation 
from  the  7  =  1.44  curve  occurs  for  Tr/Tq  «  0.  33.  This  indication  of 
condensation  occurs  at  the  spatial  location  of  X/D  *  1.  0,  as  shown  in 
Fig.  40,  beyond  which  elevated  CO2  temperatures  may  be  expected. 
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Now,  using  7  =  1.44,  the  CO2  expansion  is  plotted  in  the  P-T  plane 
along  with  its  vapor  pressure  data,  and  saturation  is  found  to  occur  at 
Tpj/T0  =  0.  57  which  is  estimated  from  Fig.  40  to  occur  at  (X/D)s  =  0.  38. 
Consequently,  CO2  expands  isentropically  fromX/D*0.4  until  reaching  an 
X/D  «  1  before  condensation  onset  occurs,  and,  assuming  a  constant 
value  for  7  and  isentropic  flow  for  the  axial  region  0.  4-^  X/D-^  1,  one 
may  calculate  the  axial  dependence  of  the  Mach  number  M  using  the 
well-known  relation: 


T/To  =11  +  [(y  -  1 ) 7 2] M 2 } 


and  the  result  is  shown  in  Fig.  54. 


(26) 


Figure  53.  C02  Rayleigh  scattering  in  the  N-T  plane,  7  *  1.44. 
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Figure  54.  Axial  variation  of  Mach  number  for  CO2. 

It  is  interesting  to  compare  the  measured  value  of  7  =  1.  44  with  the 
value  of  the  specific  heat  ratio  of  CO2  if  equilibrium  flow  occurred. 

Since  CO2  has  four  vibrational  modes,  one  of  which  is  twofold  degenerate, 
and  these  modes,  particularly  the  deformation  mode,  are  partially  ex¬ 
cited  at  284  K,  one  finds  7  of  CO2  for  this  reservoir  temperature  to  be 
1.  295.  Figure  55  shows  the  variation  with  temperature  of  7  of  CO2 
where,  for  this  calculation,  it  has  been  assumed  that  the  vibrational 
modes  are  in  equilibrium  with  themselves  and  the  translational  and  ro¬ 
tational  degrees  of  freedom.  It  is  seen  that  7  =  1.4  is  the  equilibrium 
value  for  the  temperature  region  110  K. 

The  effects  of  nonideality  of  CO2  on  the  equilibrium  specific  heat 
ratio  can  be  obtained  using  the  equation  of  state  of  the  gas.  For  this 
purpose,  the  axial  profile  of  the  static  pressure  ratio  for  the  CO2  ex¬ 
pansion  was  calculated  using  the  results  of  Fig.  53,  and  using  the 
corresponding  values  of  temperature  and  the  virial  coefficient  data  for 
the  Lennard- Jones  (6-12)  potential  (Ref.  15),  one  finds  non-ideal  gas 
corrections  to  7  for  the  CO2  expansion  whose  results  are  shown  in 
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Fig.  40  to  be  on  the  order  of  0.  5  percent;  the  calculated  values  are  shown 
in  Fig.  55.  It  should  be  noted  that  the  ideal  gas  values  of  7  for  CO2  for 
T  >  200  K  were  obtained  from  Ref.  15  and  the  values  for  T  <  200  K  were 
calculated  using  a  simple  harmonic  oscillator  model  for  each  CO2  vibra¬ 
tional  mode. 


Since  the  measured  value  of  7  exceeds  the  equilibrium  value  for 
X/D  <  1.0,  as  shown  by  Figs.  53  and  55,  it  is  reasonable  to  expect  that 
vibrational  freezing  has  occurred.  Further,  since  a  constant  7  =  1.44 
describes  the  expansion  from  X/D  ^  0.  4  to  condensation  onset,  the  vib¬ 
rational  freezing  must  have  occurred  prior  to  X/D  *  0.  4.  To  demon¬ 
strate  analytically  that  this  inference  is  reasonable  would  be  reassuring. 
For  such  a  demonstration,  AEV(12)  is  defined  as  the  vibrational  energy 
of  CO2  which  must  be  transferred  to  the  translational  and  rotational 
modes  between  axial  positions  (X/D)i  and  (X/D)2  if  the  vibrational 
modes  are  to  remain  in  equilibrium.  Since  the  vibrational  energy  is 
concentrated  in  the  deformation  modes  of  CO2  for  the  temperatures  of 
interest,  the  calculation  of  AEV(12)  can  neglect  the  symmetric  and  anti¬ 
symmetric  valence  modes  of  vibration.  If  -Jiu  is  the  quantum  energy  of 
the  deformation  mode,  AEv(12)//fiw  can  be  interpreted  by  classical  col¬ 
lision  mechanics  as  the  probability  for  transfer  of  a  vibrational  quantum 
between  points  (X/D)j  and  (X/D)2,  subject  to  the  restriction  that 
AEv(12)Mu  «  1.  The  initial  point  (X/D)^  is  chosen  to  be  the  reservoir. 
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and  (X/D>2  is  chosen  to  be  0.  3.  The  vibrational  energy  for  each  posi¬ 
tion  is  calculated  using  the  simple  harmonic  oscillator  model  of  the 
deformation  mode,  and  it  is  found  that,  between  the  reservoir  and 
X/D  =  0.  3, 


AEv(12)/Kw  =  2  x  HT2  =  pv(12) 


Next,  one  must  calculate  the  probability  pv(  12)  of  transfer  of  a 
vibrational  quantum  to  occur  between  the  axial  points  (X/D)i  and  (X/D)2. 
If  Zt(12)  is  the  average  total  number  of  collisions  which  a  molecule 
experiences  between  points  (X/D)i  and  (X/D>2  and  Ziq  is  the  number  of 
collisions  required  for  transfer  of  one  vibrational  quantum  from  vibra¬ 
tional  levels  1  to  level  0,  the  calculated  probability  of  transfer  pvc(12) 
is 


P vc(  12)  *  Zt(12)/Z10 


(27) 


If  the  molecule  is  described  insofar  as  elastic  collisions  are  concerned 
as  a  hard  sphere  of  collision  cross  section  (a^),  one  can  easily  show 
that 


(X/D)2  !/(  _]) 
Zt(12)  =  4N0a2Dv/7^y/  d(X/D)'MU  +  [(y  -  l)/2]M2i~  r 

(X/D), 


(28) 


where  all  symbols  have  the  usual  meaning.  One  sees  that  M  must  be 
known  as  a  function  of  X/D,  but  since  this  is  essentially  the  problem 
being  investigated,  for  order  of  magnitude  purposes  the  axial  variation 
of  M  for  a  7  =  1.4  expansion  is  to  be  used.  By  using  the  sonic  orifice 
density  data  for  isentropic  N2  expansions  (Ref.  9),  shown  in  Fig.  56, 
the  axial  variation  of  the  Mach  number  was  calculated  assuming  7  =  1.4, 
and  Fig.  57  shows  the  results.  With  these  results,  the  integrand  (Im) 
and  the  integral  of  Eq.  (28),  defined  as  Z^( 12 ),  were  determined  and 
are  shown  in  Fig.  58.  By  using  the  appropriate  values  of  N0,  7,  cr^, 
and  D,  it  was  found  that  Z^(0,  0.  3)  was  on  the  order  of  4000  collisions. 


65 


AEDC-TR-75-37 


X/D 

Figure  58.  Axial  variation  of  collision  number  integrand,  Im  ,  and 
accumulative  collision  number  function,  ZT(12), 

7  =  1.4(N2). 


To  determine  Ziq ,  the  data  of  Cottrell  and  McCoubrey  (Ref.  16) 
were  used,  as  well  as  the  Landau-Teller  vibrational  transfer  theory 
(Ref.  16),  to  justify  the  use  of  an  average  value  of  Ziq  of  10^  over  the 
axial  distance  range  of  interest.  Therefore, 


pvc(12)  «  4  X  10’2 


and  one  sees  that  pv(12)/pvc(12)  ~  1.  Consequently,  it  is  expected  that 
vibrational  equilibrium  may  be  maintained  for  the  axial  region 
0  <  X/D-^  0.  3.  It  is  most  reasonable  to  believe  that  vibrational  freez¬ 
ing  does  not  occur  at  the  orifice  as  the  flow  develops,  but  rather  equi¬ 
librium  or  near -equilibrium  is  maintained  to  X/D  «  0.  2  to  0.  3.  Ex¬ 
tending  point  2  farther  along  the  flow  direction  does  not  substantially 
increase  the  value  of  Z^(X/D  =  0,  X/D);  however,  the  average  value 
of  ZiQj  which  varies  as  exp  [-T"l/3],  increases  significantly,  there¬ 
by  decreasing  pvc(lzj.  Furthermore,  AEV(12)  increases  as  T  as 
point  2  decreases,  so  the  ratio  py(12)/pvc(12)  which  equals 
[AEv(12)/!iu]  /  [Zt(12)/Zio]  increases  significantly.  Obviously,  the 
meaning  is  that  the  inherent  probability  of  transfer,  as  denoted  by 
Py c ( 1 2 ) ,  is  inadequate  to  maintain  equilibrium.  Further  discussion  of 
CO2  expansions  and  condensation  must  await  additional  analysis. 
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From  the  sonic  orifice  expansion  results  for  N2  (Figs.  30  through 
33),  CO  (Fig.  36),  and  O2  (Fig.  3  7),  it  is  seen  that  the  Ashkenas- 
Sherman  prediction  for  7  =  1.4  is  in  excellent  agreement  with  the  ex¬ 
perimental  measurements  for  the  unsaturated  flow -field  regions,  and 
the  onset  of  condensation  as  manifested  by  an  increase  in  Tr  due  to 
the  recombination  heat  release  is  obvious.  Since  the  condensation  is 
the  subject  of  future  reports,  no  further  discussion  is  presented  re¬ 
garding  this  phenomenon.  However,  it  is  noted  that  the  rotational 
Raman  spectra  exhibited  a  Boltzmann  number  density  distribution  for 
the  rotational  energy  levels  for  both  the  uncondensed  and  condensed 
regions  of  the  flow  field,  thereby  demonstrating  that  local  thermal 
equilibrium  existed  for  the  rotational  mode  of  motion. 

As  Figs.  41  through  50  show,  the  conical  nozzle  expansion  results 
for  N2  density  and  rotational  temperature  along  the  centerline  are  in 
excellent  agreement  with  the  MOCS  predictions  as  the  flow  is  unsatu¬ 
rated.  Once  again,  condensation  onset  is  evidenced  by  an  increase  in 
relative  to  the  MOCS  prediction,  and  Boltzmann  distributions  were 
observed  to  adequately  characterize  the  rotational  mode  of  N2.  It  is  to 
be  noted  in  all  the  axial  profile  results  that  the  expansion  character¬ 
istics  from  the  nozzle  lip  intersect  the  centerline  at  approximately 
Xt/Df-  =  25.  Figure  59  shows  the  two-dimensional  axial-radial  locus  of 
the  expansion  characteristic,  and  obviously  the  gas  exterior  to  the  so- 
called  "testing  region"  of  the  nozzle  expands  faster  than  that  interior  to 
the  region,  thereby  resulting  in  an  earlier  saturation  of  the  gas  in  the 
exterior  region.  Consequently,  the  spatial  locus  of  saturation  is  three 
dimensional  in  nature  and  of  cylindrical  symmetry  for  the  nozzle  flow 
studied.  Rayleigh  scattering  measurements  (Ref.  10)  support  this 
assertion  and  the  previously  mentioned  inference  that  excessively  high 
off-axis  values  of  Tp  are  possibly  a  manifestation  of  condensation  even 
though  no  condensation  is  observed  along  the  centerline.  Finally,  the 
radial  profiles  of  both  the  measured  density  and  temperature  (Figs.  47 
through  50)  show  that  significant  deviations  exist  between  the  measured 
and  MOC  predicted  values  of  density  and  temperature  for  off-axis 
points  interior  to  the  testing  region.  Using  these  results,  the  experi¬ 
mentally  determined  ratio  of  the  unit  length  Reynolds  number  (Rex) 
and  Mach  number  (M)  as  a  function  of  radial  position  for  X^/Dt  =  12.47 
and  19.  79  was  found  and  is  shown  in  Fig.  60.  The  application  of  simi¬ 
lar  measurements  and  flow  field  calibrations  for  larger  flow -field  simu¬ 
lation  facilities  is  evident. 
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r/Dt 

Figure  60.  Ratio  of  the  unit  length  Reynolds  number  (Rex)  and  Mach 
number  (M)  as  a  function  of  radial  position. 
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5.0  SUMMARY  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 


5.1  REVIEW  OF  MEASUREMENTS 

A  summary  of  the  experimental  flow -field  measurements  is  shown 
in  Table  4.  The  measured  values  of  the  rotational  temperature  ranged 
from  300  to  11  K,  and  the  range  of  gas  number  density  at  which  the 
temperatures  were  measured  was  from  2.  54  x  10^  to  2.  31  x  10^ 
molecules/ cc.  At  these  respective  number  densities,  the  experimental 
uncertainty  in  Tr  was  ±2  and  ±20  percent.  Accurate  measurements  of 
the  gas  number  density  were  performed  over  the  number  density  range 
of  from  5.  7  x  10^8  to  6.0  x  10*®  molecules / cc,  corresponding  to  approxi¬ 
mately  190  to  2  torr,  respectively,  at  300  K. 


Table  4.  Summary  of  Measurements 


These  measurements  are  the  first  application  at  AEDC  of  Raman 
scattering  diagnostics  of  supersonic  and  hypersonic  flow-field  tempera¬ 
tures  and  number  density,  and  the  first  application  of  Raman  diagnos¬ 
tics  to  a  two-phase  flow  field.  The  quality  of  these  results  is  indica¬ 
tive  of  the  flow  diagnostics  capability  of  the  Raman  technique. 
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5.2  SYSTEMS  IMPROVEMENTS 

Optical  system  improvements  can  readily  increase  measurement 
accuracy  and  decrease  the  lowest  density  at  which  temperature  mea¬ 
surements  can  be  made.  Figures  61a  and  b  show  two  modifications 
which  can  be  made  to  the  present  optical  system.  Figure  61a  sche¬ 
matically  shows  a  concave  front  surface  mirror  replacing  the  laser 
dump  system  and  refocusing  the  expanding  laser  beam  back  to  the 
original  focal  volume;  i.  e.  ,  a  multiple-pass  laser  beam  system.  The 
spherical  mirror  must  have  a  focal  length  half  that  of  the  laser  focusing 
lens,  and  it  must  be  placed  a  distance  twice  its  focal  length  from  the 
scattering  or  focal  volume.  Obviously  the  problem  here  will  be  the  pre¬ 
cise  alignment  of  the  concave  mirror  and  will  require  precision  trans¬ 
lation  of  the  mirror  in  three  dimensions  and  tilting  about  the  optical 
axis.  Figure  61a  also  shows  a  concave  front  surface  mirror  replacing 
the  viewing  dump  and  placed  in  the  same  manner  as  the  laser  beam 
concave  mirror.  Thus  a  bilateral  collection  optics  system  is  created, 
and  precise  alignment  of  the  mirror  is  also  required  here.  The  beam 
double  pass  and  bilateral  collection  optics  can  increase  the  scattered 
radiation  detected  by  a  factor  of  4  (a  conservative  estimate),  and  by 
using  a  lower  f/number  collection  optics  system  than  the  one  used  in 
these  experiments,  an  additional  factor  of  4  improvement  can  be  real¬ 
ized.  Hence,  by  using  the  improvements  shown  in  Fig.  61,  scattered 
radiation  detected  can  be  increased  by  a  factor  of  16,  which  would  de¬ 
crease  the  imprecision  of  the  measurements  of  these  experiments  by  a 
factor  of  4.  Furthermore,  while  maintaining  the  precision  of  the  pres¬ 
ent  measurements  the  lower  density  limit  could  be  decreased  by  a  fac¬ 
tor  of  16  to  1.  76  x  1015  cm’l  From  another  viewpoint,  the  measure¬ 
ments  reported  herein  could  have  been  made  four  times  faster  while 
maintaining  the  same  precision. 

Figure  61b  shows  a  more  difficult  to  achieve  optical  system  improve¬ 
ment.  By  removing  the  laser  output  mirror  and  placing  it  beyond  the 
flow  field,  as  shown,  the  observed  scattering  volume  is  inside  the  laser 
cavity.  This  arrangement  requires  precise  alignment  and  excellent  sta¬ 
bility  of  the  mirror,  but  offers  the  possibility  of  a  two  order  of  magni¬ 
tude  increase  in  detected  scattered  radiation  as  shown  by  Hickman  and 
Liang  (Ref.  17).  Another  possible  problem  for  this  arrangement  is  the 
effect  of  flow -field  gradients  on  the  cavity  stability.  However,  the 
arrangement  of  Fig.  61  offers  a  possible  increase  of  a  factor  of  400  in 
the  detected  radiation  which  would  decrease  the  measurement  impre¬ 
cision  of  these  experiments  by  a  factor  of  20  and  decrease  the  lower 
density  limit  by  a  factor  of  400  to  7.  03  x  10 cm'^, 
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can  be  determined.  However,  for  a  mixture  of  such  species  as  N2  and 
CO  or  CO2,  very  high  resolution  will  be  necessary  to  separate  the  spec¬ 
tral  lines.  Since  the  centrifugal  distortion  constant  of  molecules  such 
as  N2,  O2,  CO,  and  CO2  is  very  small,  the  separation  of  the  rotational 
lines  of  a  specie  in  terms  of  wave  number  is  practically  constant. 
Therefore,  a  Fabry -Perot  interferometer  may  be  used  as  the  spectro¬ 
scopic  element,  as  shown  in  Fig.  61,  to  pass  all  the  rotational  lines 
of  a  given  specie  and  reject  the  Rayleigh  scattered  radiation  and  the 
rotational  lines  of  the  other  molecular  species  (Ref.  18).  The  total 
intensity  of  the  passed  lines  can  give  the  specie  number  density,  and 
a  small,  moderate  resolution,  single  spectrometer  can  spectrally  ana¬ 
lyze  the  output  of  the  interferometer  to  give  the  rotational  temperature 
of  the  specie. 

Laboratory  measurements  of  rotational  temperatures  and  densities 
of  single  species  at  elevated  temperatures  (>300  K)  will  be  conducted, 
and  the  effects  of  both  continuum  and  discrete  background  radiation  in 
the  spectral  region  of  the  pure  rotational  lines  on  the  accuracy  of  the 
temperature  and  density  measurements  will  be  evaluated.  These  experi¬ 
ments  will  be  extremely  important  in  determining  the  feasibility  of 
Raman  scattering  diagnostics  for  radiating  flow  fields  and  the  low- 
density  limits  for  high-temperature  measurements.  Flow -field  mea¬ 
surements  for  reservoir  temperatures  up  to  1200  K  will  also  be  made. 
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APPENDIX  A 

DETAILS  OF  TEMPERATURE  DETERMINATION 

Rotational  temperatures  are  determined  using  Eq.  (9)  and  mea¬ 
sured  rotational  line  relative  intensity  values.  If  we  set 


C 

J(J  +  1) 


b  - 


B0ho 

k  To 


and 


y  =  fn 


r  r  G  1  '  T' 

lJ  1  R 


v  -  4B  (J  +  3/2) 


)  4  |~3  (J  +  D(J  +  2)1  t 

gJ  [_2  (2 J  +  3)  J  ' 


Then  Eq.  (9)  becomes 


y  =  a  +  bx 


(29) 

(30) 

(31) 


(32) 


(33) 


For  a  least-squares  fit  of  a  straight  line  to  the  data,  a  and  b  must 
be  chosen  in  such  a  way  that  the  weighted  sum  of  the  squares  of  devi¬ 
ations  of  y  from  the  desired  straight  line  of  best  fit  is  a  minimum 
(Ref.  19). 

SvK  -  yK  -  (a  +  bxK)  (34) 


Let 


Wj  a 


2 

y  1 


v; 

K  yK 


(35) 


The  W’s  are  the  weighing  factors,  and  cr  denotes  standard  deviation.  If 
0-2  =  1,  then 


W 


K 


(36) 
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Therefore,  one  must  minimize 


K  /B  Y 

K='  \%  / 


K?1  WKSyK 


(37) 


3y£  =  yl  +  +  b24  “  “a^K  +  ^ab*K  ~  ^bx^y^ 


(38) 


KJ,  wka>'k  = 


K  K 

max  max 

£  +  ^  Wv  + 


K=  1 


K=  1 


max 

b2  2  W vxl 

1  K  K 


K  K 

max  max 


”  2a  K=i  “  bxK^  “  2b  KJj  ^KxKy K 


(39) 


(r  v,i ) 


SET 


2a  2  W*  -  2  2  Wvyv  +  2b  S  W„xv  =  0 

K=  1  K  K=  1  K  *  K=1  K  K 


(40) 


a_ 

(9b 


fe^5yi) 


SET 


=  2b  2  Wvxd  +  2a  2  Wvxv  -2  2  W^x^y*-  = 
K=  1  K  K  K=1  K  K  K=  1  R  K  R 


(41) 


K  K 

max  max 

2  Wi/yir  2  W 
K=1  K  K  K= 1 


K  K 

^  max  max 

K  XK  **  KJj  WKXK  KJj  WKxKyK 


K  K 

max  max 

2  W  v  2  Wvx£ 
K=  1  K  K=1  k  K 


-  (Sbkxic) 


(42) 


max  max 


b  = 


2  Wv  2  WvxyYy  —  S  W^x^  2  ffryv 

K=  1  K  K=1  K  K  K  K=  1  K  K  K=  1  K  K 


K  K 

max  max 

2  W„  2  Wv 
K= 1  K  K=1  K 


/«„„  v 

■b  -  U?1  *kXk/ 


(43) 


Let 


yK  = 
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From  Beer's 


K'g1). 


Then 


and  results  in 


W, 


=  fractional  standard  deviation 

(44) 

(Xi 

to 

(45) 

■  [('fir" 

(46) 

[WF 

(47) 

^ _ 

II 

(48) 

Let  Sy  be  defined  as  the  standard  deviation  of  y  about  the  line, 
From  Beer's 

IK, 

2,  SyK 


L  ^  m  a  x  \ 


(49) 


For 


Td  = 


—  h  c  B  k 


HcB/k  ,, 
ji  R  =  — — -  db 


(50) 

(51) 


Therefore, 


h<:H/k  c 

77“  bb 


(52) 
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As  stated  previously,  Eq.  (9)  is  solved  by  iterative  machine  calculation 
using  Eqs.  (30),  (31),  (32),  (42),  (43),  and  (48).  The  iteration  process 
continues  to  calculate  a  temperature  until  consecutive  calculated  values 
are  within  0.  01  percent. 

A  temperature  defined  by  the  maximum  number  of  spectral  lines  is 
first  calculated.  Then  a  temperature  value  for  Kmax  -  1  lines  is  calcu¬ 
lated.  This  process  is  repeated  i  times  through  Kmax  -  i  =  2,  i  =  0,  1, 

2,  .  .  .  .  However,  no  more  than  12  iterations  for  each  (Kmax  -  i)  num¬ 
ber  of  lines  is  allowed. 

The  standard  deviation  of  the  temperature  (Str)  is  computed  for 
each  (Kmax  -  i)  number  of  lines.  This  value  is  computed  using  Eqs.  (49), 
(52),  and  (57). 


80 


AEDC-TR-75-37 


For  a  given  flow  condition  and  flow -field  position,  more  than  one 
spectral  scan  may  be  taken.  Therefore,  it  is  necessary  to  "average" 
these  data  by  some  process.  A  portion  of  PROGRAM  RATROT  calcu¬ 
lates  an  average  Tr/Tq,  Tr/Tq,  and  an  average  standard  deviation 
(Str)  for  a  given  number  of  spectral  lines  used  in  calculating  a  Tr  in 
a  fixed  set  of  data.  The  Tr/Tq  and  Str  values  are  weighted  according 
to  their  own  precision. 


by  definition  of  a  weighted  average  (Ref.  19). 


(58) 


(59) 


(60) 


s  =  the  total  number  of  data  points  to  be  used,  usually  the  number  of 
spectral  scans  for  a  given  condition 


(61) 


by  definition  of  a  standard  deviation  of  an  average  (Ref.  16), 


(62) 


For  a  given  number  of  scans  (s)  for  a  particular  flow  condition  and 
position,  PROGRAM  RATROT  has  computed  (Tr/Tq)^  and  (<r»j*  )^.  f°r 
each  number  of  spectral  lines  used.  To  avoid  the  decision  as  to  which 
particular  (Tr/Tq)k  is  the  best  to  use,  PROGRAM  RATROT  performs 
another  weighted  average  calculation  to  determine  Tr  and  ^Tr* 


(63) 
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(64) 


(65) 
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APPENDIX  B 

OPTIMUM  ILLUMINATION  GEOMETRY  FOR  RAMAN  SCATTERING 

MEASUREMENTS 


For  the  determination  of  the  optimum  illumination  geometry  for  the 
laser-Raman  technique  for  flow  fields,  the  general  treatment  of  this  sub¬ 
ject  by  Schwiesow  (Ref.  6)  will  be  used. 

Assume  the  laser  beam  to  be  plane -parallel  of  wave  number 
~vQ  =  1/A0  and  diameter  Dl.  A  diffraction  limited  lens  of  focal  length 
Fl  focuses  the  laser  beam  in  the  flow  field.  An  axially  symmetric 
irradiance  distribution  will  exist  at  the  focus,  and  approximately  80  per¬ 
cent  of  the  total  power  in  the  incident  beam  will  be  contained  in  a  cylin¬ 
der  of  diameter 


d' 


=  (2  F  )(F,  I),  ) 


and  length 

V  ■  (U'V0)(FL'DL)2 


(66) 


(67) 


The  relative  intensity  of  the  laser  beam  at  the  focus  is 


I 


rel 

o 


I 

o 


-2n2 

^oDL 


4 


2 


(68) 


The  relative  number  of  scattered  photons  collected  per  unit  time  is 


R 


S 


V0' 


(69) 


This  can  be  rewritten  as 


Fl/Dl 

R°  “■"H—  FT)  V0' 


(70) 


The  sample  volume  (cylinder  of  diameter  d'  and  length  l')  is 
imaged  onto  the  defining  slit  of  the  spectrometer  by  the  collecting  optics. 
If  the  collection  optics  have  linear  magnification  (M),  the  sample  volume 
at  the  slit  will  have  dimensions  (d  =  Md',  JL  =  Ml')  with  the  inherent 

p 

acceptance  angle  of  the  spectrometer  enlarged  from  ^  to  =  M  ^  at 
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the  scattering  volume.  For  a  slit  of  width  w  and  height  h,  the  volume 
of  the  imaged  cylinder  is  approximately  that  of  a  rectangular  prism  d 
deep,  the  lesser  of  d  and  w  wide,  and  the  lesser  of  i  and  h  high. 

For  cases  in  which  h/w  >  i/d  over  the  possible  range  of  Fl/Dl 
values  (the  narrow  slit  case),  the  function  Rs  is  considered  over 
three  regions: 


(1) 

d  <  w, 

i  <  h 

(2) 

w  <  d. 

i  <  h 

(3) 

w  <  d. 

h  <  i 

(71) 

(72) 

(73) 


From  Eqs.  (71),  (72),  and  (73),  it  is  observed  that  the  maximum 
value  of  Rs  occurs  at  the  boundary  between  region  2  and  region  3  where 
i  =  h  =  Mi'.  At  this  point,  the  optimum  Fl/Dl  is 


1/2 


For  a  wide  slit  case,  i/d  >  h/w, 

(1)  d  <  w,  i  <  h 

(2)  d  <  w,  h  <  i 

(3)  w  <  d,  h  <  i 


the  three  regions  become 


(74) 


Regions  1  and  3  of  the  wide  slit  case  are  identical  to  regions  1 
and  3  of  the  narrow  slit  case.  In  region  2, 


R 


s2 


N. 


4\i2hn 


(75) 
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which  is  independent  of  Fl/Dl,.  Hence,  the  optimum  values  of  Fl/Dj^ 
may  be  anywhere  within  the  range  defined  by  the  boundary  of  region  2 
with  region  1  where  h  =  i  =  Mi'  and  the  boundary  of  region  2  with 
region  3  where  w  =  d  =  Md'. 


(76) 


With  a  truncated  Gaussian  laser  beam 
points,  2  changes  to  2.4  and  14  changes  to 


and  Dl  defined  by  the  1  / 
16  in  Eqs.  (74)  and  (76). 


For  the  laboratory  system. 


178.6 


and 


194.4 


for  nominal  200  Mm  by  5  mm  slits.  The  diameter  of  the  scattering  vol¬ 
ume  was  220.  5  m m.  For  the  4  by  10  Research  Chamber, 

fl 

—  =  25.5 
°L 


and 


17.4  < 


<  97.2 


for  nominal  200  Mm  by  5  mm  slits.  The  diameter  of  the  scattering  vol¬ 
ume  was  31.5  Mm  with  a  length  of  5.  35  cm. 
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NOMENCLATURE 


a 


Nucleon  number  of  the  atomic  constituent 


B0 

b 


C 

Cp 


CRR 

c 

D 


Rotational  constant  of  the  ground  vibrational  level 

Defined  by  Eq.  (31)  as  the  slope  of  the  rotational 
Boltzmann  plot 

A  constant  defined  by  Eq.  (5) 

A  calibration  factor  for  number  density  measurement  and 
defined  by  Eq.  (12) 

A  constant  as  defined  in  Eq.  (2) 

Speed  of  light  in  vacuum 

Sonic  orifice  diameter 


Dt 

d',d 


FL 

gj 

h. 


JA>  JB 


I 


J' 


Diameter  of  laser  beam  before  focusing 
Throat  diameter  of  conical  nozzle 

Diameter  of  the  scattering  volume  and  the  diameter  of 
the  scattering  volume  as  imaged  on  the  spectrometer 
slit,  respectively 

Focal  length  of  the  laser  beam  focusing  lens 

Nuclear  spin  statistics  weighting  factor 

Planck's  constant  or  spectrometer  slit  height;  Ti  =  h/27r 

Nuclear  spins  for  species  A  and  B,  respectively 

Absolute  and  relative  intensity  (photons/sec)  of  a  pure 
rotational  Raman  line 


% 

lo. 


lo 


I 


T' 


Collision  number  integrand 

Laser  beam  intensity  (photons/sec-cm^)  before  and  after 
focus,  respectively 

Absolute  total  and  relative  total  intensity  of  the  rotational 
Raman  Stokes'  band 


J  Rotational  quantum  number 

K,Kmax  A  summation  index  (1,2,3,  .  .  .)  indicating  the  rotational 

spectral  line  number  and  the  maximum  number  of  spectral 
lines  used  in  a  rotational  temperature  determination, 
respectively 
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k 

O e 

M 

MOC 
N,  Ng,  Nc 
No 
o 
P 

P  P 
ro>  r 

p 

Pv<12> 


qr 


ci 

r 


Rex 

Hs 


r/Dt 

SjR 

Sj 

str.sTr 


Sb 


Boltzmann's  constant 

Length  of  the  scattering  volume  and  the  length  of  the 
scattering  volume  as  imaged  on  the  spectrometer  slit, 
respectively 

Magnification  of  the  collection  optics  system  or  Mach 
number 

Method  of  characteristics 

Molecular  species  number  density,  cm'^ 

Reservoir  gas  number  density,  cm"2 
Ortho  modification 
Electronic  parity  eigenvalue 

Reservoir  pressure  and  gas  static  pressure,  respectively 
Para  modification 

Probability  of  transfer  of  a  vibrational  quantum  between 
(XfD)1  and  (X/D)2 

Molecular  rotational  partition  function 

Quantum  mechanical  and  classical  molecular  rotational 
partition  function,  respectively 

Unit  length  Reynolds  number 

Relative  number  of  scattered  photons  collected  per  unit 
time 

Nondimensional  radial  distance  from  flow  centerline 

Fractional  standard  deviation  of  the  relative  intensity  of  a 
spectral  line 

Rotational  Raman  line  strength  factor 

Standard  deviation  of  the  determined  rotational  tempera¬ 
ture  as  defined  by  Eq.  (52)  and  average  standard  devi¬ 
ation  of  the  determined  rotational  temperature  as  defined 
by  Eq.  (64),  respectively 

Standard  deviation  of  the  slope  of  the  straight  line  fitted  to 
the  data 


S 


y 


Standard  deviation  of  the  data  about  the  fitted  straight  line 
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Tj  Tc 

T°  _ 
TR>TR 

Tr(J) 

T 
A  s 

V 

Wj 

WK 

WK 

w 

X/D,  Xt/Dt 

(X/D)s 

x 

y^k 

yj 

yk 

z10 

ZT(12) 

Zt(12) 

7 

AEV(12) 

0R 

"o 


Molecular  species  temperature,  K 
Reservoir  gas  temperature,  K 

Rotational  temperature  and  average  rotational  tempera¬ 
ture  as  defined  by  Eq.  (63),  respectively 

Sensitivity  factor  of  optics,  spectrometer,  and  detector 

Static  temperature  at  saturation 

Scattering  volume 

A  weighting  factor  as  defined  by  Eq.  (62) 

Rotational  line  intensity  weighting  factor 
A  weighting  factor  as  defined  by  Eq.  (65) 

Spectrometer  slit  width 

Nondimens  ion  al  axial  distance  from  sonic  orifice  and 
conical  nozzle  throat,  respectively 

Axial  location  of  saturation 

J(J  +  1) 

Defined  by  Eq.  (32) 

Defined  by  Eq.  (60) 

Defined  by  Eq.  (58) 

Number  of  collisions  required  for  transfer  of  one  vib¬ 
rational  quantum  from  vibrational  level  1  to  level  0 

Average  total  number  of  collisions  which  a  molecule 
experiences  between  (X/D)^  and  (X/D>2 

Accumulative  collision  number  function 

Ratio  of  specific  heats 

Vibrational  energy  of  CO2  transferred  to  the  translational 
and  rotational  modes  between  (X/D)i  and  (X/D) 2 

Characteristic  rotational  temperature 

Wavelength  of  rotational  Raman  line 

Laser  line  frequency 

Laser  line  wave  number 
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P>  P o 

<J 

CT(lKrel) 

aj 

aR 

^RRj 

(^tr)k 

T 

u 


Gas  density  and  reservoir  gas  density,  respectively 
Hard  sphere  collision  cross  section 

Standard  deviation  of  the  relative  intensity  of  a  rotational 
Raman  line 

Defined  by  Eq.  (59) 

Total  rotational  Raman  scattering  cross  section 

a 

Rotational  Raman  line  cross  section,  cm^/ster 

Standard  deviation  of  the  average  Tr  determined  using 
K  spectral  lines 

Time 

Collection  optics  solid  angle  and  spectrometer  solid  angle, 
respectively 

Vibrational  frequency  of  CO2 
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